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1.  Abstract 

We  have  theoretically  and  experimentally  investigated  the  dynamics  of  frequency  comb 
generation  in  silicon-based  microresonators.  In  this  system,  the  parametric  gain  associated  with 
the  nonlinear  process  of  four-wave  mixing  leads  to  parametric  oscillation  and  generation  of 
signal/idler  sidemodes.  This  report  describes  our  results  that  further  develop  the  parametric 
comb  technology.  We  theoretically  investigate  comb  generation  and  modelocking  dynamics,  and 
the  connection  to  concepts  of  synchronization  and  self-organization.  We  also  investigated  mode- 
crossings  in  microresonators  and  tunable  coupling,  and  their  impact  on  comb  formation.  We 
studied  comb  generation  without  an  external  pump  laser.  In  addition,  we  explored  comb 
generation  at  other  pump  wavelengths  using  dispersion  engineering.  Finally,  we  investigated 
comb  generation  via  coherent  supercontinuum  generation  in  silicon  nitride  waveguides. 

Our  research  effort  illustrates  that  the  silicon-based  platform  offers  a  high  potential  as  a 
platform  for  robust,  integrated,  chip-scale  comb  source.  Since  the  dispersion  of  the  system  can  be 
readily  tailored  through  waveguide  cross-section  engineering,  the  platform  offers  flexibility  in 
operating  wavelength  and  bandwidth.  The  developments  made  here  are  a  significant  step  towards 
the  development  of  chip-based  frequency  combs  that  can  operate  in  real-world  sensing  and 
spectroscopy,  for  optical  clocks,  and  for  frequency  metrology  applications. 
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2.  Introduction 

Over  the  past  decade,  there  has  been  tremendous  development  in  the  development  of 
techniques  for  producing  optical  frequency  combs  [1],  with  applications  including  spectroscopy, 
optical  clocks,  arbitrary  waveform  generation,  frequency  metrology,  and  astronomical 
spectrograph  calibration  [2,3,4].  Traditionally,  modelocked  solid-state  and  fiber  lasers  have  been 
used  for  optical  frequency  comb  generation.  For  full  self-stabilization  of  frequency  combs,  it  is 
desirable  to  have  a  broadband,  octave-spanning  spectrum  for  stabilizing  the  carrier-envelope 
offset  frequency /ceo,  which  is  one  of  the  two  degrees  of  freedom.  Using  an  f-2f  interferometer, 
the  low-frequency  component  of  the  comb  is  doubled  and  heterodyned  with  the  high-frequency 
component,  allowing  for  extraction  of  the  /ceo,  which  can  be  locked  to  a  microwave  frequency 
standard  [5].  The  other  degree  of  freedom  is  the  comb  spacing,  or  the  repetition  frequency /rep, 
which  can  be  directly  detected  and  referenced  to  a  frequency  standard,  such  that  the  mth  comb 
line  has  frequency  fm  =/ceo  +  «?//?• 

However,  a  novel  approach  to  producing  frequency  combs  has  emerged  in  the  past  10  years 
through  the  use  of  parametric  nonlinear  interactions  in  microresonators  [6,7].  In  this  approach,  a 
single  CW  laser  beam  coupled  to  a  microresonators  can  produce  stabilized,  octave- spanning 
combs  through  highly  cascaded  four-wave  mixing  (FWM)  oscillation  and  soliton  modelocking. 
The  tight  transverse  confinement  of  the  light  within  these  resonators  enables  operation  at  modest 
CW  pump  powers  and  powerful  dispersion  engineering  that  allows  the  for  the  creation  of 
anomalous  dispersion  over  a  wide  range  of  wavelength  ranges  necessary  to  produce 
ultrabroadband  combs.  The  highly  compact  size  of  the  this  comb-generating  system  ensures  that 
it  will  operate  at  low  powers  to  produce  combs  with  spacings  ranging  from  10  GHz  to  1  THz 
while  being  highly  integrated  and  robust  and  an  enabling  technology  for  many  spectroscopic 
applications  and  devices.  In  our  research,  we  utilize  silicon  nitride  microresonators  for 
parametric  frequency  comb  generation  [8,9,10,11].  Unlike  many  other  high-<2  resonator  designs, 
the  resonator  and  the  coupling  waveguide  are  monolithically  integrated.  Thus,  the  entire  on-chip 
configuration  of  CMOS -compatible  microresonators  can  provide  robust  operation  in  ambient 
conditions  yielding  a  truly  monolithic,  robust,  and  highly  compact  frequency  comb  source. 
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3.  Modeling  Parametric  Frequency  Comb  Generation 

We  use  the  Lugiato-Lefever  (LL)  model  [12,13,14,15,16]  to  simulate  the  full  temporal 
dynamics  of  comb  generation  dynamics  in  microresonators.  We  extend  this  model  to  include 
higher-order  dispersion  and  self-steepening  to  enable  simulations  of  comb  spectra  that  span  an 
octave.  The  LL  model  has  demonstrated  excellent  agreement  with  experimental  results  when 
used  to  solve  for  a  steady-state  solution  of  the  spectrum  [14].  In  our  analysis  we  assume  the 
mode  frequencies  are  equally  spaced,  using  the  free-spectral  range  (FSR)  of  the  pumped  mode, 
and  allow  the  phase  of  these  modes  to  define  the  dispersion-induced  frequency  offset. 
Additionally,  since  we  treat  the  field  both  in  the  time  and  frequency  domains,  additional  terms 
such  as  higher-order  dispersion,  self-steepening,  Raman  scattering,  two-  and  three-photon 
absorption,  and  free-carrier  effects  can  easily  be  included  in  a  manner  similar  to  supercontinuum 
simulations  [17].  Using  this  model,  we  identify  distinct  and  necessary  stages,  including  one  that 
is  chaotic  and  unstabilized,  that  yields  a  path  to  stable,  single-pulse  modelocking. 

We  express  the  LL  equation  with  contributions  from  higher-order  dispersion  and  self- 
steepening  as,  where  E  is  the  resonator  field,  Em  is  the  pump  field  at  frequency  coo,  t  and  r  are  the 
slow  and  fast  times,  respectively,  describing  the  slow  evolution  of  the  cavity  field  and  temporal 
field  within  the  cavity  at  a  given  time  t,  y  is  the  nonlinear  parameter,  L  is  cavity  length,  tR  is  the 
round-trip  time,  So  is  the  cavity  detuning,  a  is  loss  per  roundtrip,  6  is  the  transmission  coefficient 
between  the  resonator  and  the  bus  waveguide,  and  /? n  is  the  n-th  order  dispersion  coefficient. 


dE(t,  t) 
dt 


+  iyL 


1  +  - 


co0  dr  j 


\ 

\E(t,r)f 


E(t,r)  +  JdEm 


(1) 


Figure  1  shows  a  route  to  achieving  stable,  single-pulse  modelocked  octave- spanning 
frequency  comb.  Parameters  are  taken  from  the  226-GHz  FSR  Si;,N4  microring  resonator  in 
Ref.  [10].  The  continuous-wave  (CW)  pump  wavelength  lp  at  1550  nm  is  positioned  in  the 
anomalous  group-velocity  dispersion  (GVD)  regime  between  the  two  zero-GVD  points 
(ZGVD’s)  of  the  waveguide  ZGVD1  and  ZGVD2  at  1000  nm  and  1760  nm,  respectively,  with  a 

I  |2 

pump  power  Pm  =|£'in|  of  1.5  W.  The  system  goes  through  four  stages,  corresponding  to  four 

different  frequency  detunings  of  the  pump  field.  Stage  I  consists  of  the  formation  of  multiple 
stable  cavity  solitons.  At  the  bare-cavity  detuning  So  =  0,  the  resonance  shifts  away  from  the 
pump  due  to  the  Kerr  effect  as  energy  builds  up  within  the  ring.  Anomalous  GVD  allows 
modulation  instability  to  convert  energy  from  the  driven  mode  to  the  primary  offset  resonances 
~23  modes  (50  nm)  away,  which  is  dictated  by  phase-matching  conditions  [18].  Further  build-up 
of  power  then  causes  cascaded  four-wave  mixing  (FWM)  of  these  MI  peaks  within  the 
anomalous  GVD  region.  Various  other  combination  of  FWM  processes  will  also  occur,  which 
results  in  the  formation  of  mini-combs.  Once  these  mini-combs  begin  to  overlap,  the  temporal 
characteristics  change  dramatically,  forming  stable  cavity  solitons  as  initial  modelocking  occurs. 
In  Stage  II,  as  the  pump  is  tuned  closer  to  the  Kerr-shifted  resonance,  the  intracavity  power 
increases,  and  the  multi-pulse  modelocking  is  lost  due  to  instability  of  the  cavity  solitons. 
Adjacent  pulses  interact  with  each  other  causing  destabilization  and  what  appears  to  be  chaotic 
behavior  [19].  Further  detuning  of  the  pump  past  the  Kerr- shifted  resonance  in  Stage  III  causes 
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t  [ns] 


the  intracavity  power  to  drop,  and  modelocking  is  again  observed.  After  the  change  in  So, 
multiple  cavity  solitons  re-form,  but  only  the  stable  ones  persist.  The  remaining  pulses  collide 
with  the  nearest  neighbors  and  annihilate,  causing  spikes  in  power  and  broad  lines  seen  in  the 
corresponding  spectral  evolution.  In  the  final  stage,  the  pump  detuning  is  increased  again.  At  the 
lower  intracavity  powers  the  remaining  cavity  solitons  destabilize  and  dissipate.  At  the  suitable 
pump  detuning,  the  system  evolves  to  a  single  pulse  within  the  microresonator  The  modelocked 
frequency  comb  is  smooth  and  spans  an  octave  across  ZGVD2  into  normal  GVD  regime,  as  has 
been  observed  experimentally  [10],  with  the  dispersive  wave  feeding  energy  into  the  longer 
wavelength  modes. 
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Figure  1.  Evolution  of  single-pulse  modelocking  within  a  226  GHz  monolithic  Si3N4  ring 
resonator.  Each  of  the  four  stages  (bottom  to  top)  represent  an  increase  in  pump  detuning  from  on- 
resonance  to  So  =  0.02,  0.04,  0.05642  at  simulated  times  r  =  25,  50,  75  ns.  The  pump  power, 

\E-m\  is  kept  constant  at  1.5  W.  Column  1  is  the  temporal  evolution,  with  power  plotted  on 


identical  linear  scales,  and  I  spanning  a  single  roundtrip  time.  Column  2  is  the  spectral  evolution, 
with  spectral  power  plotted  on  the  same  100  dB-scale.  Column  3  and  4  show  the  temporal  and 
spectral  power,  respectively,  at  the  end  of  Stages  I  -  IV. 


Next,  we  investigate  the  effects  of  GVD,  higher-order  dispersion,  and  self-steepening  on 
sideband  generation  from  FWM  [20]  by  analyzing  the  coupled  mode  equations  associated  with 
the  field  E(t ,  r)  =  A0  +  A+  +  A_  [21,22,23],  where  Ao  is  the  pump  field  and  A+  and  A_  represent  the 

symmetrically  detuned  sidemodes.  For  our  analysis,  we  assume  that  the  amplitude  of  the 
sidebands  are  much  smaller  than  Ao,  in  which  case  the  coupled  equations  are  given  as, 


tR  —  =  ~(a+i60)A0 +iyL\A0f  A0  +  2iyLA*0A+A_+JdEin, 
dt 


(2) 
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Power  in  mode  [dB] 


(3) 


tR^  =  -  a  +  iAK+iL  Y  ^Qn  +2iL— y\A0\2  A++iyL(l~— )ao2A 
9t  [  „=15....  n[  «o  J  "o/ 


tR^  =  -  a  +  iAx-iL  Y  OlQ"  -2iL—y\A0\2  A  +iyL\l  +— }a20A*+, 
9t  [  „=X5,...  Ul  M0  J  OJ0 

where  Q  is  the  sideband  detuning,  and  the  phase  mismatch  is  given  by, 


Ak=S„-L  J  Aor-2yL\A„ 

«=2, 4,...  n' 

The  resulting  solutions  for  A+  and  A.  are, 


(4) 


(5) 


(6) 


where  A(£2)  =  -a  +  ^yL(l-£22/n>o)|A0|2J  -(A k)2  is  the  gain  coefficient.  This  gain 

coefficient  A  and  phase  mismatch  A k  depend  on  the  even-orders  of  dispersion,  while  the  odd- 
order  dispersion  terms  add  a  phase  contribution  to  the  sidebands  which  is  opposite  in  sign  with 
respect  to  each  other.  While  the  odd-order  terms  do  not  contribute  directly  to  the  amplitude 
growth,  these  phase  terms  are  relevant  for  pulse  shaping  and  mode-locking.  The  inclusion  of  the 
self- steepening  term  does  not  affect  A k,  but  modifies  the  nonlinear  term  in  the  gain  and 
contributes  an  additional  phase  term. 

We  investigate  how  the  comb  bandwidth  can  be  controlled  using  two  different  waveguide 
cross  sections  of  690x1900  nm  and  910x1800  nm.  Previously,  the  high  stress  in  silicon-nitride 
films  prevented  deposition  of  high  quality  nitride  films,  which  in  turn  limited  the  magnitude  of 
the  anomalous  GVD  that  could  be  reached.  However,  recent  advances  in  fabrication  techniques 
have  allowed  for  deposition  of  thick,  high  quality  films  [24],  enabling  higher  anomalous  GVD 
values  and  access  to  a  wider  parameter  space  for  waveguide  dispersion.  The  GVD  for  the  S^lSU 
waveguide  is  shown  in  Fig.  2(a)  which  was  calculated  with  a  finite-element  mode-solver.  For  the 
comb  simulations,  the  microring  has  a  100-pm  radius  and  is  pumped  at  1550  nm  with  650  mW  in 
the  coupling  waveguide.  Figure  2(b)  shows  the  generated  comb  spectra  for  the  two  different 
dimensions.  While  the  region  of  anomalous  GVD  for  the  690x1900  nm  cross  section  spans  only 
550  nm,  the  low  anomalous  GVD  at  the  pump  wavelength  allows  the  comb  to  span  a  bandwidth 
far  exceeding  the  anomalous  GVD  region.  In  contrast,  for  the  910x1800  nm  cross  section,  for 
which  the  anomalous  GVD  region  spans  1250  nm,  the  large  anomalous  GVD  at  the  pump 
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wavelength  results  in  narrowband  comb  generation.  We  find  that  the  overall  comb  bandwidth 
that  is  ultimately  produced  is  dictated  in  part  by  the  frequency  detuning  between  the  pump  and 
the  first  generated  sideband.  Considering  the  phase  mismatch  equation  with  only  the  GVD 
contribution,  the  detuning  is  inversely  proportional  to  the  square  root  of  [h.  This  implies  that  for 
larger  (smaller)  [h  values  the  smaller  (larger)  the  frequency  detuning  results  in  a  narrower 
(broader)  bandwidth.  The  bandwidth  dependence  on  GVD  is  consistent  with  the  predictions 
based  on  steady-state  analysis  by  Coen  and  Erkintalo  [25].  Figure  2(c)  shows  the  temporal 
profile  of  the  simulated  combs  and  indicates  single-pulse  modelocking  and  cavity  soliton 
formation,  which  is  indicative  of  a  stabilized  comb  [11,26].  The  narrower  bandwidth  allows  for  a 
higher  power  per  comb  line  allowing  for  more  efficient  power  conversion  from  the  pump  to  the 
comb  lines,  which  is  critical  for  an  energy-efficient  multiple  wavelength  source  for  WDM 
applications. 
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Figure  2.  (a)  Simulated  GVD  for  Si3N4  waveguides  with  cross  sections  690^1900  nm  (blue) 
and  910x1800  nm  (green),  (b)  Simulated  comb  spectra  and  (c)  temporal  plot  for  260-GHz  FSR 
Si3N4  microring  resonators  with  cross  sections  690x1900  nm  (top)  and  910x1800  nm  (bottom). 
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Figure  3.  Measured  comb  spectra  in  Si3N4  microrings  with  cross  sections  (a)  690x1900  nm 
and  (b)  910x1800  nm  (green). 


We  investigate  these  predictions  experimentally  by  pumping  at  a  microring  resonance  near 
1560  nm  using  an  amplified  single-frequency  laser.  To  cover  the  entire  spectral  measurement 
range,  we  utilize  two  optical  spectrum  analyzers  operating  from  900  -  1200  nm  and  from  1200  - 
2400  nm.  Figures  3(a)  and  3(b)  show  the  measured  spectra  from  the  690x1900  nm  and 
910x1800  nm  cross-section  microring,  respectively.  The  690x1900  nm  microring  [Fig.  3(a)] 
generates  a  158-THz  bandwidth  comb,  which  represents  the  broadest  comb  generated  to  date  in  a 
silicon  nitride  platform.  The  additional  spectral  lines  between  2  -  2.4  pm  [shaded  region  in  Fig. 
3(a)]  correspond  to  second-order  diffraction  of  the  comb  lines  generated  in  the  1-1.2  pm  range 
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in  the  grating-based  spectrometer.  In  contrast,  the  comb  generated  in  the  690x1900  nm  microring 
[Fig.  4(a)]  spans  a  significantly  narrower  bandwidth  of  28  THz.  The  bandwidths  of  the  measured 
spectra  are  in  good  agreement  with  our  theoretical  predictions  [Fig.  2(b)]. 

Next,  we  investigate  the  effects  of  higher-order  dispersion  by  modeling  comb  generation  at 
three  different  pump  wavelengths,  1550  nm,  1400  nm,  and  1200  nm.  Figure  4(a)  shows  the  GVD 
and  fourth-order  dispersion  (FOD)  for  a  690x1900  nm  cross  section,  and  Fig.  4(b)  shows  the 
simulated  comb  spectra  for  each  pump  wavelength.  In  each  case,  a  stabilized,  singlepulsed 
modelocked  comb  is  achieved.  In  comparing  the  comb  spectra  for  the  1550-  and  1400-nm 
pumps,  the  comb  bandwidth  is  significantly  broader  for  the  1400-nm  pump  even  though  the 
magnitude  of  GVD  is  smaller  for  the  1550-nm  pump.  This  is  due  to  the  FOD  contribution  to 
phase  matching,  which  becomes  significant  for  FWM  over  broader  bandwidths.  The  broadest 
comb  is  generated  with  the  1200-nm  pump,  where  the  GVD  and  FOD  values  are  the  lowest. 
While  the  GVD  at  1200  nm  (/?2  =  -16  ps2/km)  is  close  to  that  at  1550  nm  (fh=  -17  ps2/km),  as  a 
result  of  the  different  FOD  values  at  1200  nm  (Jh  =  2.6  x  10'4  ps4/km)  and  1500  nm  (Jm=  1.8  x 
10'3ps4/km),  the  30-dB  bandwidth  with  a  1200-nm  pump  is  1.8  times  broader  than  that  with  a 
1550  nm  pump.  Our  results  show  that,  while  low  anomalous  GVD  is  a  necessary  requirement  for 
broadband  combs,  small  FOD  is  also  critical  to  further  extend  the  bandwidth. 
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Figure  4.  (a)  Simulated  GVD  and  FOD  for  690x1900  nm  Si3N4  waveguide,  (b)  Simulated 
comb  spectra  for  pump  wavelengths  of  1550  nm,  1400  nm,  and  1200  nm. 


In  addition,  we  observe  dispersive  wave  (DW)  formation  at  the  edges  of  the  comb  spectra 
[Fig.  4(b)].  Recent  studies  have  investigated  the  role  of  third-order  dispersion  (TOD)  and  FOD 
on  DW  formation  [25,27].  The  wavelength  for  DW  formation  can  be  predicted  by  phase¬ 
matching  conditions  between  the  cavity  soliton  and  the  DW  wave  across  a  ZGVD  point  [28]. 
While  TOD  and  FOD  contribute  to  DW  generation,  the  group  velocity  of  the  soliton  and  the 
DW,  along  with  the  higher-order  dispersion  terms  must  be  taken  into  account  to  accurately 
predict  the  spectral  position  for  DW  generation. 
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4.  Self-Organization  in  Soliton  Modelocked  Parametric  Frequency 
Combs 

A  large  collection  of  coupled  oscillators  with  slightly  different  natural  frequencies  can 
undergo  a  transition  to  a  phase-locked  state  with  identical  frequencies.  This  phenomenon  appears 
in  many  systems  spanning  biology,  chemistry,  neuroscience,  and  physics  [29,30].  Examples 
include  power  grid  networks,  neural  networks,  chemical  oscillators,  and  arrays  of  Josephson 
junctions  and  semiconductor  lasers  [31,32,33,34,35].  Self-organization  in  such  systems  has  been 
modeled  by  the  Kuramoto  model,  which  describes  the  time-evolution  of  the  phase  </>p(t)  of  an 
oscillator  p  as  an  interaction  between  its  natural  frequency  cop  and  its  coupling  to  the  phases  of 

all  the  other  oscillators.  The  governing  equations  are  <pp  -  cop  +  K ^  sin(^m  ,  where  k  is  the 

coupling  strength  [30].  This  model  can  be  recast  in  an  order-parameter  formulation,  where  an 

average  phase  and  a  coherence  R(t)  are  defined  via  R(t)t"v  =— V  e4"'  .  Then  the  Kuramoto 

model  becomes  <j>p  =  (op  +  KR(t)sinUp-<pp^ .  Viewed  this  way,  (f>p  is  no  longer  coupled  to  every 

individual  oscillator's  phase,  but  only  to  the  average  phase  if/.  Moreover  the  effective  strength  of 
the  coupling  is  proportional  to  the  coherence  R.  This  proportionality  between  coupling  and 
coherence  creates  a  positive  feedback  which,  for  a  sufficiently  large  k,  gives  rise  to  an  abrupt 
transition  in  which  a  macroscopic  fraction  of  the  oscillators'  frequencies  spontaneously 
synchronize. 

In  optics  an  alternative  form  of  phase  locking  can  occur  in  lasers  and  parametric  oscillators 
between  a  large  collection  of  cavity  modes  with  nearly  equidistant  frequency  separations.  In 
these  systems  the  nearest  neighbor  spacing  across  the  modes  varies  due  to  dispersion  within  the 
cavity.  In  the  presence  of  nonlinearity  within  the  cavity,  the  system  can  spontaneously  mode- 
lock  such  that  the  frequency  spacings  between  the  oscillating  modes  become  identical  (Fig. 
lc,d).  Although  Kerr-based  parametric  frequency  combs  have  been  suggested  as  the  most 
fundamental  example  of  self-organization  in  nonlinear  optics  [36,37],  no  direct  connection  has 
been  made  to  the  concepts  of  synchronization  and  self-organization.  We  find  that  phase 
equations  derived  from  the  Lugiato-Lefever  equation  (LLE)  display  self  organization  features 
akin  to  the  Kuramoto  model,  including  the  existence  of  meaningful  order  parameters  and 
coherence-coupling  feedback  that  describes  the  soliton  formation  process  [38].  These  equations 
predict  that  the  pump  phase  is  offset  from  the  rest  of  the  modal  phase  profile  in  the  soliton  state. 
Additionally  our  analysis  predicts  that  phase  anti-symmetrization,  where  the  phase  profile 
becomes  anti-symmetric  about  the  pump  phase,  occurs  before  phase  synchronization  and  soliton 
formation  can  occur. 

Using  the  LL  equation  and  employing  a  slowly  varying  envelope  approximation  for  the  total 
intra-cavity  field  and  normalized  dispersion  coefficients  =(ljivg/L\k  vgfik ,  we  derive  the 
following  general  dynamical  phase  equation,  with  time  dependency  of  the  phases  made  implicit: 
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(8) 


&  &  Nl 2 

^=yp2+fpJ-r  2 

^  l,m,n=-N/2 

where  T  =yL/tR  and  A'"ip  =AlAmAn/ Ap  .  This  equation  has  functional  similarities  to  the 

Kuramoto  model  wherein  each  optical  mode  can  be  considered  an  individual  oscillator.  The 
spread  in  natural  frequencies  of  the  oscillators  is  represented  by  the  second-  and  third-order 
dispersion  terms  in  the  right-hand  side,  while  the  nonlinear  term  gives  rise  to  the  coupling  among 
oscillators  in  the  last  term  in  the  right-hand  side.  In  the  absence  of  a  strong  pump  mode,  this 
equation  has  no  stables  solutions.  We  further  simplify  this  equation  by  keeping  only  terms  that 
correspond  to  the  pump-degenerate  (PD)  FWM  processes,  where  two  pump  photons  are 
annihilated  to  created  a  photon  pair,  and  the  pump-nondegenerate  (PND)  FWM  processes,  where 
one  pump  photon  and  one  comb  photon  are  annihilated  to  create  two  photons.  From  this 
approximation,  we  obtain  a  pair  of  equations  that  correspond  to  the  phase  average 
<f)p  =  (<pp  +  <p_p) /2  and  difference  dp  =(<pp  -(p_p)/2p  for  pairs  of  modes  symmetric  about  the  pump 

mode.  We  term  the  following  the  parametric  synchronization  equations  (PSE): 

$P  =jP2  -2rr72A^cos[2(0o -^p)]-TriA^NR(t)cos(^  -jp)cos[p(0p -d0)\  (9) 


0P  =^rP2  ~  2TriAcN  R(t)sm(<p0 -(pp)sin[P(dp  -fl0)], 

j  p 


(10) 


where  //  is  the  ratio  between  a  pump  and  comb  mode,  Ac  =A0/rj  is  the  amplitude  of  the  comb 

modes,  <f>o  is  the  pump  phase  which  is  fixed  and  N  is  the  total  number  of  comb  modes  in  the 
system.  The  phase  average  equation  describes  the  symmetric  behavior  of  the  system  while  the 
phase-difference  equation  describes  its  anti-symmetric  behaviorThe  pump-nondegenerate 
processes  have  both  symmetric  and  anti-symmetric  contributions  and  appear  as  the  last  terms  in 


the  phase-average 
8  v'N2 


and  phase  difference  equations.  R(t)  — 


N 


'N2  im(em-e.) 

e 

/  m=l 


and 


90(t )  = — -  y  m6m  represent  the  order  parameters  of  the  system.  Here  60  is  the  normalized 

m=A 

average  phase  difference.  It  serves  the  same  role  as  the  average  phase  t//  in  the  Kuramoto  model. 


One  of  the  key  predictions  of  the  PSE  is  that  the  system  evolves  into  a  state  in  which  the 
pump  phase  is  offset  from  the  rest  of  the  phase  profile.  This  offset  arises  from  the  cosinusoidal 
dependence  of  the  PD  term;  in  order  for  this  term  to  act  as  a  restoring  force  on  the  phase  average 
as  in  the  Kuramoto  model,  it  must  have  a  sinusoidal  dependence  on  the  phase  averages.  We 
compare  in  detail  the  synchronized  spectral  phase  profiles  in  Fig.  5.  Both  systems  stabilize  to  a 
broadband  phase-locked  state  with  an  offset  of  the  pump  phase  from  the  rest  of  the  phase  profile. 
Due  to  the  factor  of  2  in  the  argument  of  the  PD  term,  this  offset  should  be  0  <  ^offset  <  tt/2  and 
centered  at  tt/4  in  order  for  the  cosine  to  have  a  significant  sine-like  contribution.  Both  the  PSE 
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and  the  LLE  predict  pump  phase  offsets  within  these  bounds.  The  LLE  system  has  a  slightly 
larger  pump  phase  offset  due  to  self-phase  and  pump-induced  cross-phase  modulation  effects 
that  were  not  accounted  for  in  the  PSE.  This  is,  to  our  knowledge,  the  first  theoretical  prediction 
and  explanation  of  the  origin  of  the  pump  phase  offset  of  the  soliton-modelocked  states  in  a 
parametric  frequency  comb.  We  confirm  in  Fig.  5(c,d  insets)  that  the  broadband  phase-locked 
state  results  in  a  solitary  pulse  in  the  time  domain.  The  exact  pulse  shape  for  the  PSE  is  not 
quantitatively  meaningful  since  all  the  modes  have  equal  amplitude  resulting  in  a  Sine -like  pulse 
without  a  CW  background. 
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Figure  5.  Three  stages  of  evolution  of  the  phase  profile  of  the  intra-cavity  field  as  predicted  by 
(a)  the  parametric  synchronization  equtions  (PSE)  and  (b)  the  Lugiato-Lefever  equations  (LLE): 
The  grey  curves  represent  the  initial  random  phase  profile.  The  red-dotted  curves  show  the  phase 
evolution  after  370  (PSE),  308  (LLE)  round  trips;  the  phase  profiles  in  both  models  illustrate  the 
anti-symmetrization  of  the  spectral  phase  due  to  the  pump  degenerate  LWM  processes.  The  black 
curves  show  the  final  phase  profile  after  3394  (both)  round  trips  in  which  the  phases  have  become 
completely  synchronized,  which  is  a  result  of  the  pump  non-degenerate  terms.  In  addition,  a  slight 
offset  of  the  phase  of  the  pump  from  the  phases  of  the  other  cavity  modes  is  observed.  (c,d)  The 
final  spectral  phase  profiles  of  the  PSE  and  LLE  systems  showing  deviations  from  a  pure  linear 
profle,  including  the  pump  phase  offsets,  (insets)  Temporal  pulse  shapes  of  the  PSE  and  LLE. 
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5.  Mode  Coupling  in  Microresonators 

The  microresonator  resonances  can  in  principle  be  precisely  calculated  using  the  dispersion 
of  the  resonating  modes  and  the  resonator  length.  However,  modal  coupling  between  different 
types  of  modes  can  significantly  alter  the  shape  and  position  of  their  resonances.  Mode  splitting 
occurs  for  strong  coupling  [39],  and  coupling  between  whole  families  of  modes  results  in 
avoided  crossings  [40,41,42,43,44,45].  Avoided  crossings  occur  in  resonators  when  two  distinct, 
frequency-degenerate  modes  are  coupled  and  form  new  eigenmodes.  The  frequencies  of  these 
eigenmodes  are  no  longer  degenerate  with  a  splitting  proportional  to  the  coupling  strength 
resulting  in  the  absence  of  a  crossing  of  the  eigenmode  resonances.  In  microresonators  this 
behavior  can  lead  to  dramatic  localized  changes  in  the  effective  dispersion  near  these  crossing 
points,  which  in  general  affects  any  parametric  interaction  that  relies  on  precise  frequency 
matching  of  different  resonances.  In  particular  it  can  play  an  important  role  in  the  formation  of 
parametric  frequency  combs  [41,42,43,44,45,46,47,48].  In  the  context  of  frequency  comb 
generation,  only  modal  interactions  between  different  families  of  spatial  modes  have  been 
considered  thus  far.  However,  in  dielectric  waveguides,  even  when  the  waveguide  is  'single 
mode',  there  are  typically  at  least  two  guided  fundamental  modes,  the  fundamental  quasi 
transverse  electric  (TE0o)  and  the  fundamental  quasi-transverse  magnetic  (TM0o)  mode,  which 
correspond  approximately  to  the  polarization  of  light  in  the  waveguide.  Since  such  a  mode 
interaction  can  even  occur  in  single-mode  waveguides,  it  is  more  universal  than  other  forms  of 
modal  interactions  (i.e.,  between  higher-order  spatial  modes).  The  physical  origin  and  strength  of 
the  modal  coupling  between  the  TEoo  and  TMoo  modes  are  based  on  different  parameters  of  the 
ring  resonator  such  as  its  radius  of  curvature,  waveguide  cross-section,  and  side-wall  angle 
[49,50].  While  in  most  cases  the  polarization  properties  follow  from  the  corresponding  straight- 
waveguide  characteristics  [51],  microresonators  with  smaller  bend  radii  and  larger  side-wall 
angles  typically  will  exhibit  greater  modal  coupling. 
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Figure  6.  (a)  Schematic  of  the  setup  used  for  observing  and  characterizing  polarization  mode 
coupling  in  microresonators  on  a  temperature  stabilized  (T)  silicon  chip.  (b),(c)  SEM  images  of 
the  microresonators  and  the  coupling  regions  with  the  bus  waveguides  for  (b)  the  1.8-mm-long 
microresonator  and  (c)  a  100-pm  radius  ring  resonator. 

Our  experimental  setup  for  investigating  polarization  mode  coupling  is  depicted  in  Fig.  6. 
We  probe  the  resonators  with  two  different  external-cavity  diode  lasers  covering  a  total  tuning 
range  between  1450  nm  and  1640  nm.  Lensed  fibers  are  used  to  couple  into  and  out  of  the  bus 
waveguide  with  inverted  tapers  for  mode-matching.  The  polarization  of  the  input  and  output  light 
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is  controlled  and  analyzed  with  standard  fiber-based  polarization  controllers  and  a  polarization 
beam  splitter,  and  the  output  power  is  monitored  with  a  sensitive  photodiode.  We  use  a 
temperature  controller  with  a  Peltier  element  on  the  chip  holder  to  stabilize  and  tune  the  Si3N4 
microring  resonators  under  investigation.  To  overcome  the  limited  precision  of  our  tunable 
lasers,  we  use  an  automated  stepped  scanning  and  fitting  routine  supplemented  by  calibrating 
each  resonance  position  with  a  high-precision  wavemeter.  We  find  that  this  method  leads  to  an 
average  precision  better  than  50  MHz. 
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Figure  7.  Optical  transmission  for  TM0o  polarized  light  and  selected  resonances  (a)  -  (e)  in 
detail  for  TM0o  (blue)  and  TE0o  (orange).  The  appearance  of  the  TE0o  resonances  for  TM0o 
polarized  input  (and  vice  versa)  and  their  behavior  indicate  an  avoided  crossing  near  1595  nm. 


We  first  investigate  polarization  mode  coupling  in  two  ST,N4  microrings  with  725x1100  nm 
and  725x900  nm  waveguide  cross-sections  (effective  refractive  index  of  1.7  for  both),  lOO.pm 
radii  and  10°  side-wall  angles  [52],  The  transmission  measurement  for  TMoo  input  light  for  the 
first  microring  (Fig.  7)  yields  sharp  and  deep  resonances  due  to  nearly  critical  coupling  between 
the  bus  waveguide  and  the  resonator.  Near  1595  nm  a  second  sharp  resonance,  which  we 
associate  with  the  TEoo  mode,  appears  on  the  right  side  of  the  main  TM0o  resonance  and  becomes 
deeper  until  both  show  the  same  extinction.  The  main  resonance  then  experiences  an  adiabatic 
crossover  and  the  secondary  resonance  (now  on  the  left  side)  slowly  disappears.  We  attribute  this 
behavior  to  an  avoided  crossing  at  1595  nm  associated  with  a  strong  modal  interaction  between 
the  TEoo  and  TM0o  modes. 

We  further  investigate  the  effects  of  both  polarization  and  higher-order  mode  crossings  on 
comb  generation  by  directly  comparing  the  comb  spectra  generated  in  microring  resonators  with 
a  measurement  of  its  mode-crossings  (Fig.  8).  The  combs  are  generated  by  strongly  pumping  a 
single  resonance  at  1540  nm.  We  investigate  a  microring  resonator  with  a  waveguide  cross- 
section  of  725x1650  nm2  (effective  refractive  index  of  1.8),  a  100-pm  radius,  and  a  10°  side-wall 
angle  [Fig.  8(a)].  We  observe  a  polarization  mode-crossing  near  1580  nm  with  a  splitting  at  the 
avoided  crossing  of  2.2  GHz,  which  greatly  exceeds  the  intrinsic  loss  rate  of  the  resonator.  This 
polarization  mode  crossing  manifests  as  a  reduction  in  the  mode  intensity  in  the  generated  comb 
spectrum.  In  addition,  there  exists  an  anomaly  in  the  FSR  at  1550  nm,  which  can  be  attributed  to 
a  mode  crossing  with  a  higher-order  spatial  mode,  and  we  observe  a  similar  corresponding 
feature  in  the  comb  spectrum.  Furthermore,  we  observe  suppressed  comb  generation  when 
pumping  near  one  of  the  higher  order  mode  crossings  which  we  attribute  to  the  large  change  in 
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the  effective  dispersion.  Since  this  effect  can  be  disadvantageous  for  applications  related  to 
frequency  comb  generation,  it  needs  to  be  taken  into  account  in  the  optimal  device  design. 
Moreover,  we  investigate  the  effects  of  mode  crossings  on  comb  generation  in  an  80-GHz  FSR 
microresonator  with  a  waveguide  cross-section  of  725x1700  nm2  (effective  refractive  index  of 
1.8),  a  1.8-mm  length,  and  a  10°  side-wall  angle  close  [Fig.  8(b)].  For  this  microresonator  our 
measurements  reveal  three  polarization  mode  crossings  and  a  number  of  higher-order  mode 
crossings  that  affect  the  comb  spectrum  generated  from  this  resonator. 


Figure  8.  Generated  comb  spectra  in  Si3N4  microresonators.  Insets  show  the  measured  FSR’s 
with  the  (upper)  blue  branch  being  for  the  TE00  mode  and  revealing  the  mode  crossings,  (a) 
Results  for  a  725x1650  nm  cross-section,  100-pm  radius  microring.  A  polarization  mode  crossing 
(pol)  occurs  near  1580  nm  with  a  corresponding  feature  (indicated  by  an  arrow)  in  the  comb 
spectrum.  A  second  anomaly  in  the  FSR’s  near  1550  nm  due  to  a  mode  crossing  with  a  higher- 
order  spatial  mode  (horn)  also  produced  a  corresponding  feature  in  the  spectrum,  (b)  Results  for  a 
725x1700  nm  cross-section  1.8-mm  length  microresonator.  A  number  of  polarization  and  higher 
order  mode  crossings  affect  the  generated  comb  spectrum. 


We  identify  different  strategies  that  could  minimize  the  disruptive  effects  of  mode-crossings 
(polarization  and  higher-order  mode)  on  frequency  comb  generation.  In  general,  mode-crossings 
are  especially  disruptive  when  they  appear  directly  at  the  pump  wavelength  or  when  they  line  up 
symmetrically  with  respect  to  the  pump  since  they  will  then  simultaneously  affect  both  the  signal 
and  idler  resonances.  Since  the  position  of  any  mode-crossing  depends  sensitively  on  the  exact 
length  of  the  microresonators,  slight  variations  in  the  design  of  the  resonator  lengths  can  be  used 
to  circumvent  both  scenarios.  In  addition,  if  the  target  application  allows  it,  increasing  the  FSR 
will  reduce  the  frequency  at  which  mode-crossings  will  occur.  Another  approach  for  reducing 
mode-crossing  disruptions  would  be  to  minimize  the  modal  interactions  altogether,  which  may 
be  achieved  by  optimizing  the  microresonator  design  and  fabrication,  including  larger  bend  radii 
or  smaller  side-wall  angles.  Furthermore,  we  have  observed  experimentally  that  higher-order 
mode-crossing  are  suppressed  for  the  narrower  and  thus  more  symmetric  cross-section  resonators 
as  compared  to  the  wider,  more  asymmetric  cross-section  resonators.  This  indicates  that  the 
aspect  ratio  can  be  used  as  an  additional  parameter  for  controlling  mode-crossing  effects. 
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6.  Tunable  Frequency  Combs  using  Dual-Microring  Resonators 

Tuning  the  coupling  between  the  bus  waveguide  and  resonator  is  challenging  since  in 
standard  passive  monolithically-integrated  structures,  the  coupling  gap  is  fixed  by  design  and 
cannot  be  changed  after  fabrication.  Demonstrations  of  tunable  comb  generation  have  included 
resonance  frequency  tunability,  which  enables  control  over  the  operating  wavelength.  Such 
tuning  has  been  achieved  during  comb  generation  via  both  thermal  and  electro-optic  means 
[6,53,54],  Coupling  gap  tunability  enables  greater  control  over  the  comb  generation  process, 
allowing  for  optimization  of  comb  efficiency  [55].  Coupling  tunability  has  been  demonstrated  in 
integrated  silicon  devices  using  a  Mach-Zehnder  interferometer  (MZI)  coupler  [56,57,58]. 
However,  the  coupling  produced  by  this  MZI  structure  inherently  has  a  sinusoidal  wavelength 
dependence,  which  can  be  detrimental  to  broadband  comb  generation.  Coupling  gap  tunability 
can  occur  in  suspended  microtoroid  and  wedge  resonator  structures,  since  the  devices  are 
operated  using  suspended  tapered  fibers  [6];  however,  this  approach  requires  highly  stable  fiber 
positioning  with  accuracy  on  the  scale  of  tens  of  nanometers  in  order  to  achieve  controllable 
tuning.  We  show  a  dual-cavity  coupled  microresonator  structure  in  which  tuning  one  microring 
resonance  frequency  induces  a  change  in  the  overall  cavity  coupling  condition,  as  evident  in  the 
transmission  extinction  ratio  [59]. 


Figure  9.  (a)  Schemetic  of  dual-cavity  coupled  microring  resonator  with  integrated  micro¬ 
heaters  (b)  Simulated  mode  anti-crossing  curve  as  a  function  of  effective  mode  index  detuning 


between  the  two  microrings.  Ring  #2  is  tuned  while  ring  #1  is  kept  constant.  Insets  above  are 
simulated  transmission  spectra  showing  varying  extinction  across  the  anti-crossing. 

Our  structure  consists  of  two  identical  microring  resonators  evanescently  coupled  to  each 
other,  with  one  ring  coupled  to  a  bus  waveguide  [Fig.  9(a)].  Due  to  the  evanescent  coupling,  the 
two  resonant  cavity  modes  hybridize  to  form  a  coupled  eigenmodes  system,  in  which  two 
superposition  eigenmodes  (“supermodes”)  are  formed:  a  symmetric  (s)  and  an  anti-symmetric 
(as)  mode.  These  two  new  eigenmodes  exhibit  modified  resonant  frequencies,  according  to: 


(11) 


in  which  ft)avgis  the  average  of  the  individual  cavity  resonance  frequencies,  A co  is  the  difference 
between  the  individual  cavity  resonances  (cavity  detuning),  and  k0)  is  the  inter-ring  temporal 
coupling  rate.  If  the  cavities  are  degenerate  (Acu  =  0),  the  supermode  resonances  are  split  apart 
from  the  isolated  cavity  resonance  frequency  by  ±k0).  When  the  inter-ring  coupling  rate  exceeds 
the  individual  cavity  decay  rate,  distinct  doublet  resonances  form  at  each  cavity  free  spectral 
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range  (FSR).  As  the  detuning  between  the  two  cavities  is  varied,  a  characteristic  anticrossing 
shape  is  formed,  as  in  Fig.  9(b).  At  zero  detuning,  light  resonantly  couples  between  the  two 
rings,  resulting  in  an  equal  distribution  of  light  in  both  cavities  [Fig.  9(b),  middle  inset],  and  thus 
equal  extinction  in  both  doublet  resonances.  However,  as  the  cavity  detuning  increases,  the 
distribution  of  light  becomes  unequal  as  the  resonant  coupling  between  cavities  becomes  less 
efficient.  Far  from  zero  detuning,  the  modes  are  less  hybridized,  and  most  of  the  light  in  each  of 
the  doublet  resonances  is  concentrated  in  only  one  ring  or  the  other  [Fig.  9(b),  left/right  insets]. 
The  change  in  detuning  also  causes  the  effective  cavity  length  to  change,  which  alters  the 
balance  between  the  round-trip  loss  and  the  bus-waveguide  coupling.  This  modifies  the  loaded- 
Q.  The  intrinsic- 2  remains  unchanged,  as  long  as  both  cavities  have  the  same  loss  rate.  Figure 
9(b)  shows  the  simulated  doublet  resonance  shape  for  three  detuning  positions  along  the  anti¬ 
crossing  curve,  showing  the  strong  asymmetry  in  extinction  away  from  zero-detuning. 
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Figure  10.  (a)  Micrograph  of  fabricated  dual-cavity  device  with  integrated  platinum  micro¬ 
heaters.  (b)  Measured  (blue  points)  mode  anti-crossing  curves  with  theoretical  curve  fit  (red  line). 
Ring  #2  is  tuned  while  ring  #1  is  kept  fixed.  The  inset  infrared  micrographs  show  spatial  light 
distribution  among  the  two  rings  for  3  different  positions  along  anti-crossing,  (c)  Experimental 
measurement  of  extinction  vs.  detuning  for  a  critically-coupled  device  and  (d)  for  an  over-coupled 
device.  Ring  #2  is  detuned  while  ring  #1  is  kept  fixed,  (e)  Experimental  measurement  of  extinction 
tuning  at  a  fixed  wavelength  by  using  heater  #1  for  compensation.  As  heater  #2  is  increased, 
heater  #1  is  decreased  in  order  to  keep  the  resonance  wavelength  fixed. 


Our  fabricated  devices  consist  of  coupled  high-Q  SbN4  microresonators  with  tunable 
extinction  ratio  via  integrated  microheaters.  We  fabricate  a  dual  coupled  microresonator  in  SbN4 
with  cross-sectional  dimensions  950x1400  nm.  Our  fabricated  devices  have  an  intrinsic  2-factor 
of  approximately  2xl06,  and  radii  of  115  pm  as  well  as  75  pm,  yielding  devices  with  FSR’s  of 
200  GHz  and  500  GHz,  respectively.  This  cross-section  yields  anomalous  GVD  (fii  =  -180 
ps2/km)  at  our  pump  wavelength  of  1560  nm  critical  to  ensure  phase  matching  for  parametric 
comb  generation  [60].  We  do  not  observe  a  decrease  in  intrinsic-Q  compared  to  equivalent 
single-ring  devices;  this  can  enable  equivalent  comb  performance  in  the  dual-ring  device.  Above 
the  waveguide  cladding,  we  fabricate  integrated  microheaters  by  sputtering  platinum  and  using  a 
lift-off  approach,  yielding  a  cross-section  of  100  nm  tall  by  6  mm  wide.  We  position  the  heaters 
1.9  mm  above  the  waveguide  to  ensure  negligible  optical  loss  while  maintaining  close  proximity 
for  efficient  heat  delivery.  These  heaters  yield  an  efficiency  of  1.35x1 0'5  RlU/mW.  We  can 
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control  the  heater  power  to  within  7  mW,  resultinxxxxxg  in  a  resonance  wavelength  accuracy  of 
7  MHz,  which  is  well  below  the  cavity  linewidth. 

A  micrograph  of  the  200  GHz  FSR  device  is  shown  in  Fig.  10(a).  Using  the  integrated 
microheaters,  we  first  tune  ring  #2  while  keeping  ring  #1  constant.  In  order  to  achieve  blue 
detuning  (An  <  0),  we  bias  heater  #1  above  room  temperature.  The  blue  data  points  in  Fig.  10(b) 
are  measured  resonance  positions  as  the  heater  on  ring  #2  is  tuned.  The  red  line  in  Fig.  10(b)  is  a 
curve  fit  based  on  Eq.  (11),  with  an  additional  parameter  included  to  account  for  thermal  cross¬ 
talk  between  the  two  rings.  In  Figs.  10(c)  and  (d),  we  plot  the  measured  extinction  as  a  function 
of  detuning  for  a  critically-coupled  device  and  an  over-coupled  device,  respectively.  Figure  10(e) 
shows  our  ability  to  compensate  for  the  overall  wavelength  shift  as  the  extinction  is  thermally 
tuned.  By  tuning  both  heaters  independently  for  each  ring,  we  are  able  to  compensate  for  this 
wavelength  shift,  and  demonstrate  13.3  dB  of  tuning  of  the  resonance  extinction,  from  0.7  dB  to 
14  dB.  Since  we  can  keep  the  resonance  in  place  as  the  extinction  is  tuned,  this  device  is  useful 
for  real-world  applications  involving  single-frequency  lasers. 
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Figure  11.  (a)  Transmission  measurement  of  500  GHz  FSR  dual-coupled  microring  resonator. 
Measured  FSR  is  shown  for  (b)  lower  wavelength  resonance  and  (c)  higher  wavelength  resonance. 
The  corresponding  loaded  Q  is  shown  for  (d)  lower  and  (e)  higher  wavelength  resonance.  The 
generated  comb  is  pumped  at  (f)  1560  nm  on  a  left  resonance  and  (g)  1540  nm  on  a  right 
resonance,  indicated  by  solid  black  vertical  lines.  The  shaded  region  indicates  locations  of  mode- 
crossings  that  cause  degradation  of  the  comb  line  power. 


The  wide  tunability  of  our  device  allows  us  to  overcome  mode-crossings  resulting  from 
strong  coupling  between  different  transverse  spatial  modes  -  a  major  challenge  for  comb 
generation.  When  mode-crossings  occur,  one  expects  the  FSR  to  deviate  strongly  from  the 
expected  value,  accompanied  by  a  reduction  in  the  Q  caused  by  the  significantly  enhanced  losses 
of  the  higher-order  modes.  In  order  to  characterize  the  mode-crossings  in  the  dual-cavity 
structure,  we  measure  the  FSR  and  loaded  Q  of  both  supermodes  of  a  500  GHz  FSR  dual-cavity 
device.  The  transmission  measurement  of  the  dual-coupled  microring  is  shown  in  Fig.  11(a)  and 
the  measured  FSR  and  the  loaded  Q  factor  are  shown  in  Figs.  1 1(b)— 1 1(e).  The  resonance 
wavelengths  are  precisely  determined  using  a  laser-based  precision  measurement  of  the 
wavelength-dependent  FSR,  which  allows  for  measurement  of  the  FSR  with  a  relative  precision 
of  10'4  [52],  The  measurement  indicates  that  the  presence  of  mode-crossings  severely  disrupts  the 
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resonance  frequency  position.  For  the  left  resonances,  we  observe  large  higher-order  mode- 
crossings  at  1570  nm,  1585  nm,  and  1610  nm,  which  result  in  significant  deviations  in  the  FSR 
and  additional  small  mode-crossings  at  1520  nm  and  1545  nm.  For  the  right  resonances,  we 
observe  large  higher-order  mode-crossings  at  1560  nm  and  1585  nm  and  a  small  mode-crossing 
at  1525  nm.  These  regions  are  indicated  by  gray  sections  in  Figs.  ll(b)-(g).  The  corresponding 
loaded  Q  characterization  shows  that  the  large  mode-crossings  are  accompanied  by  a  significant 
reduction  in  the  Q  caused  by  the  significantly  enhanced  losses  of  the  higher-order  modes.  As  a 
result,  pumping  at  these  wavelengths  hinders  the  generation  of  a  stable  comb  due  to  insufficient 
power  enhancement  in  the  cavity.  However,  as  shown  in  Fig.  11(f)  and  11(g),  by  ensuring  that 
the  spectral  position  of  these  mode  crossings  is  far  detuned  from  the  pump  at  1560  nm  and  1540 
nm  for  the  left  resonances  and  right  resonances,  respectively,  we  are  able  to  generate  a  stable, 
low-  comb  with  1.5  W  of  pump  power.  The  generated  spectra  show  some  asymmetry,  largely 
due  to  the  presence  of  the  large  mode-crossings.  In  addition  to  altering  the  spectral  position  of 
the  resonance,  these  mode-crossings  reduce  the  cavity  enhancement  and  subsequently  restrict  the 
comb  bandwidth.  We  demonstrate  that  this  device  can  be  used  to  avoid  mode-crossings  by 
dynamically  tuning  the  position  of  a  mode-crossing.  The  ability  to  dynamically  control  the 
position  of  mode-crossings  is  vital  particularly  for  many  applications  that  cannot  rely  on  the 
tuning  of  the  pump  source.  Beginning  with  an  as-fabricated  200  GHz  FSR  dual-cavity  device 
with  100  mW  of  heating  on  ring  #1  [Fig.  12,  blue  curve],  we  observe  a  mode-crossing  at  1548 
nm,  indicated  by  a  local  decrease  in  extinction  for  a  single  FSR  resonance.  When  we  sweep  the 
resonance  of  ring  #2  by  0.4  nm  (heater  power  from  0  to  40  mW),  we  observe  that  the  mode¬ 
crossing  shifts  across  the  spectrum  by  over  3  nm  (2  FSR’s),  which  is  almost  an  order  of 
magnitude  larger  than  the  ring  resonance  shift.  Therefore,  we  can  tune  the  mode-crossing  with 
high-efficiency.  This  high  tuning  efficiency  is  due  to  the  anti-crossing  behavior  of  the  dual  cavity 
structure.  This  novel  degree  of  freedom  can  be  used  for  comb  generation  optimization. 
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Figure  12.  Demonstration  of  mode-crossing  tenability  using  a  200  GHz  FSR  dual-cavity 
device.  Ring  #1  is  heated  with  100  mW  and  ring  #2  is  swept  with  heater  power  indicated  in  the 
legend.  A  mode-crossing  (circled)  is  tuned  across  two  FSR’s  while  the  overall  position  of  the 
desired  resonances  shifts  only  0.4  nm. 
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Figure  13.  Comb  generation  dynamics  for  500  GHz  FSR  microresonator  as  pump  wavelength 
is  tuned  into  resonance  (top  to  bottom).  Plot  shows  (a)  optical  spectra  (b)  RF  spectra,  and  (c)  eye 
diagram  of  single  comb  line  that  is  filtered  and  modulated  with  a  10  Gb/s  PRBS. 


Furthermore,  for  the  generated  comb  to  be  used  as  a  multiple  wavelength  WDM  source,  low 
RF  amplitude  noise  is  required.  We  investigate  the  comb  generation  dynamics  and  the  properties, 
in  particular,  to  determine  the  stability  of  a  comb  line.  We  pump  the  resonator  at  1560  nm  for 
comb  generation.  To  monitor  the  comb  generation  dynamics,  we  measure  the  optical  and  RF 
spectra,  as  well  as  the  eye  diagram  through  modulation  of  a  single  comb  line.  For  the  eye 
diagram,  we  use  a  1  nm  tunable  filter  to  pick-off  a  single  comb  line  which  is  modulated  with  a 
231-1  non-retum-to-zero  pseudo-random  bit  sequence  (NRZ  PRBS)  at  10  Gb/s  and  sent  to  a  high¬ 
speed  sampling  oscilloscope  for  characterization.  Figure  13(a)-(c)  shows  the  measured  optical 
spectra,  RF  spectra,  and  the  eye  diagram,  respectively,  as  the  pump  wavelength  is  tuned  into 
resonance  (top  to  bottom).  While  the  optical  spectra  show  no  significant  changes,  the  RF  spectra 
show  that  the  comb  undergoes  a  transition  to  a  low  RF  amplitude  noise  state,  similar  to  behavior 
observed  previously  [11].  Furthermore  the  eye  diagram  illustrates  the  behavior  of  the  comb  at  a 
single  resonance.  In  the  high-noise  state  (top,  middle),  the  eye  shows  poor  signal-to-noise  and 
significant  distortion,  which  can  be  attributed  to  the  fast  intensity  fluctuations  of  multiple  comb 
lines  within  the  single  resonance  through  various  FWM  processes.  As  the  pump  is  tuned  further, 
the  signal-to-noise  further  degrades  to  the  point  where  the  eye  is  completely  closed,  which 
corresponds  to  comb  instability  from  chaotic-like  behavior  [16].  Once  the  comb  transitions  to  the 
low  phase  noise  state,  the  eye  diagram  shows  good  signal-to-noise.  Our  results  indicate  that, 
even  with  mode-crossings  in  close  proximity  in  wavelength,  it  is  possible  to  generate  a  stable, 
low-noise  comb  suitable  as  a  multiple  wavelength  source  for  WDM  applications. 
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7.  Microresonator-Based  Comb  Generation  Without  an  External 
Laser  Source 

Frequency  comb  generation  in  microresonator  platforms  requires  tuning  an  external,  single¬ 
frequency,  continuous-wave  (cw)  laser  into  a  resonance  of  the  microresonator.  As  the  pump 
power  is  coupled  into  the  microresonator,  thermal  effects  shift  the  resonance  to  higher 
wavelengths,  creating  a  soft  thermal  lock  between  the  cavity  resonance  and  the  pump  laser  [61]. 
When  the  intracavity  power  exceeds  the  threshold  for  parametric  oscillation,  cascaded  four-wave 
mixing  and  higher-order  four-wave  mixing  processes  occur,  resulting  in  comb  generation. 
Fluctuations  in  the  frequency  or  power  of  the  pump  laser  can  detune  the  pump  out  of  the  cavity 
resonance,  which  disrupts  comb  generation.  These  limitations  can  be  overcome  using  a  drop-port 
of  the  microresonator  as  shown  by  Peccianti  et  al.  [62,63]  with  a  hydex  microring  imbedded  in  a 
fiber  cavity  pumped  with  an  erbium-doped  fiber  amplifier  (EDFA).  Through  filter-driven  four- 
wave  mixing,  they  demonstrate  a  high-repetition-rate  pulse  source  with  a  60-nm  bandwidth  and  a 
200-GHz  repetition  rate.  In  this  drop-port  configuration  the  microresonator  acts  as  a  filter, 
seeding  the  EDFA  with  wavelengths  that  correspond  to  resonances  of  the  microresonator,  similar 
to  the  recent  work  done  by  Cholan,  et  al.  [64]  for  realization  of  a  multiple  wavelength  source. 


FPC  EDFA  Bandpass  Filter 

Figure  14.  Experimental  setup  for  comb  generation.  The  system  is  based  on  a  dual-cavity 
design  consisting  of  the  silicon-nitride  microresonator  and  an  external  fiber  cavity. 

We  explore  a  dual-cavity  architecture  where  a  single  bus  waveguide  (through-port),  which 
forms  part  of  the  external  fiber  cavity,  is  coupled  to  the  microresonator  (Fig.  14).  Amplified 
spontaneous  emission  (ASE)  from  an  EDFA  is  coupled  into  a  silicon-nitride  microresonator.  A 
fiber  polarization  controller  (FPC)  allows  for  adjustment  of  the  polarization  of  light  coupled  to 
the  microresonator,  which  is  a  critical  issue  for  cw  lasing  of  the  external  fiber  cavity  and  comb 
generation.  A  polarizer  is  placed  at  the  output  to  select  quasi-TE  polarization.  The  dispersion  of 
the  microresonator  waveguide  is  engineered  to  optimize  parametric  FWM  gain  and  oscillation 
for  the  TE  modes.  Due  to  the  higher  density  of  photonic  states  in  the  microresonator  as  compared 
to  the  fiber  cavity,  preferential  emission  occurs  at  the  microresonator  modes  resulting  in  lasing 
and  parametric  comb  generation  defined  by  the  microresonator.  Our  through-port  configuration 
is  effectively  analogous  to  microresonator-based  comb  generation  using  a  cw  pump  laser. 
Additionally,  since  only  a  single  bus  waveguide  is  coupled  to  the  microresonator,  the  system 
operates  with  reduced  coupling  losses  and  allows  for  higher  power  efficiency.  Through  the  use 
of  this  dual-cavity  design,  we  achieve  broadband  comb  generation  spanning  more  than  880  nm. 
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Key  advantages  of  this  scheme  are  that  the  system  requires  only  a  narrow-band  optical  amplifier, 
as  opposed  to  a  stabilized  single-frequency  laser,  and  that  it  eliminates  the  need  to  pump  at  and 
tune  to  a  resonance  wavelength  which  eliminates  disruptions  to  comb  generation  due  to  pump 
frequency  fluctuations  which  can  shift  the  pump  out  of  resonance. 

We  apply  the  dual-cavity  design  to  produce  a  broad  comb  spectra  using  a  230-GHz  free 
spectral  range  (FSR)  microresonator  with  a  waveguide  cross  section  of  725x1600  nm  to  provide 
a  broad  region  of  anomalous  dispersion  which  yields  a  comb  spectrum  spanning  880  nm  (90 
THz)  with  over  390  comb  lines  when  pumped  with  2.17  W  of  EDFA  power  [Fig.  15(a)].  For  a 
microresonator  with  an  80-GHz  FSR  and  725x1700  nm  cross  section  a  bandwidth  of  730  nm  (78 
THz)  [Fig.  15(b)]  is  achieved  with  2  W  of  EDFA  power. 


Wavelength  (nm) 

Figure  15.  Parametric  frequency  comb  spectra  generated  from  silicon-nitride  microresonators 
with  FSR’s  of  (a)  230-GHz  and  (b)  80-GHz. 

In  order  to  understand  the  spectral  properties  of  our  dual-comb  system,  we  perform 
simultaneous  RF  and  optical  spectral  measurements  of  the  dual-cavity  output.  We  utilize  the 
dual-cavity  configuration  with  the  bandpass  filter  removed  from  the  external  fiber  cavity.  By 
adjusting  the  polarization,  we  allow  the  external  fiber  cavity  to  lase  without  interaction  with  the 
resonator  [see  Fig.  16(a)].  The  output  from  the  bus  waveguide  is  detected  with  a  12.5-GHz 
photodiode,  and  the  resulting  RF  signal  is  amplified  with  a  0. 1  -  500-MHz  low-noise  amplifier 
and  measured  using  an  RF  spectrum  analyzer.  The  detected  RF  beatnote  is  5  MHz  [Fig.  16(c)], 
which  corresponds  to  beating  of  the  external  cavity  modes  since  the  1-GHz  line  width  of  the 
microresonator  mode  supports  multiple  modes  of  the  external  cavity.  The  inset  in  Fig.  20(c) 
further  illustrates  this  multi-mode  behavior  as  the  RF  noise  is  present  out  to  2  GHz.  Next,  the 
polarization  in  the  bus  waveguide  is  rotated  to  quasi-TE  for  comb  generation  and  the  entire  comb 
output  is  sent  to  the  photodiode  and  RF  amplifier  for  characterization.  The  measured  optical 
power  at  the  photodiode  is  equal  to  that  for  the  case  of  lasing  in  the  fiber  cavity.  When  a  state  of 
steady  comb  generation  is  reached  [Fig.  16(b)],  despite  the  multi-mode  pump,  a  steady  state  of 
comb  generation  can  be  reached  in  which  the  RF  amplitude  of  the  5  MHz  beatnote  is  reduced  by 
20  dB,  and  the  beatnote  linewidth  broadens  significantly  [Fig.  16(d)].  In  addition,  Fig.  16(d) 
inset  shows  220  MHz  (3-dB  bandwidth)  of  RF  noise  during  comb  generation,  indicating  that  the 
number  of  oscillating  external  cavity  modes  is  greatly  reduced  when  comb  generation  occurs. 
We  believe  the  system  can  be  further  improved  to  allow  for  single-mode  operation  by  reducing 
the  external  cavity  length.  Additionally,  we  expect  that  complete  stabilization  of  the  comb  can  be 
achieved  through  FSR  control  of  both  the  microresonator  and  external  fiber  cavity. 
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Microresonator  FSR  stabilization  can  be  implemented  through  pump  power  or  temperature 
control,  while  stabilization  of  the  external  fiber  cavity  FSR  can  be  performed  using  a  piezo-based 
delay  arm. 
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Frequency  (MHz) 


Figure  16.  Optical  output  spectra  without  (a)  and  with  (b)  coupling  to  the  microresonator  and 
(c)  and  (d)  are  the  corresponding  RF  spectra,  respectively.  The  insets  in  (c)  and  (d)  show  the  RF 
noise  from  DC  to  2  GHz. 
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8.  Frequency  Comb  Generation  Near  Visible  Wavelengths 

Comb  generation  near  the  visible  wavelength  range  is  of  particular  interest  since  there  are 
several  atomic  absorption  lines  such  as  rubidium  and  cesium,  to  which  frequency  combs  can  be 
locked,  providing  stable  frequency  references  for  frequency  metrology  and  optical  clockwork.  In 
addition,  the  800  to  900  nm  range  is  of  interest  for  biological  applications  since  it  is  a  suitable 
compromise  between  increased  scattering  at  shorter  wavelengths  and  increased  water  absorption 
at  longer  wavelengths.  The  high  coherence  of  frequency  comb  lines  also  allows  for  deeper  tissue 
penetration  and  higher  resolution  imaging,  particularly  for  optical  coherence  tomography  [65]. 
Previous  attempts  at  broadband  comb  generation  near  the  visible  wavelength  range  have  resulted 
in  combs  that  are  either  narrowband  or  fall  below  the  detection  noise  floor  well  before 
approaching  visible  wavelengths.  Comb  generation  closer  to  the  visible  wavelength  range  has 
been  previously  demonstrated  by  pumping  at  1  pm  [66],  but  the  shortest  wavelength  to  which 
this  comb  could  generate  was  950  nm.  These  combs  fail  to  extend  to  shorter  wavelengths  since 
the  GVD  of  these  microresonator  designs  cannot  adequately  compensate  for  the  stronger  material 
dispersion  at  shorter  wavelengths.  Anomalous  GVD  is  required  for  the  FWM  gain  that  produces 
parametric  oscillation,  but  the  increasingly  strong  normal  GVD  of  the  material  as  the  wavelength 
approaches  the  absorption  band  edge  is  fundamental  to  most  optically  transparent  materials. 

First,  we  demonstrate  broadband  comb  generation  with  a  1-pm  pump.  We  engineer  the 
dispersion  of  a  ST,N4  ring  resonator  to  overcome  the  stronger  material  dispersion  at  lower 
wavelengths,  and  we  simulate  the  comb  formation  dynamics  to  ensure  broadband  comb 
generation.  We  use  anomalous  waveguide  dispersion  due  to  the  high  confinement  of  the  optical 
mode  to  adequately  compensate  for  the  normal  material  GVD  of  S13N4.  While  previous 
demonstrations  have  also  used  this  approach,  we  further  optimize  the  dispersion  of  the  ring 
resonators  by  simulating  the  comb  formation  dynamics  for  various  cross  sectional  shapes.  We 
simulate  the  comb  spectrum  [shown  in  Fig.  17(a)]  using  the  modified  LL  model  [16].  The 
microresonator  has  a  radius  of  46  pm  and  is  pumped  at  1060  nm  with  1.5  W  in  the  coupling 
waveguide.  Our  simulations  indicate  that  the  comb  output  exhibits  stable,  single-pulse 
modelocking  behavior.  Modelocking  gives  precisely  uniform  comb  spacing,  which  is  desirable 
for  many  of  the  metrology  and  clocking  applications  of  frequency  combs.  Experimentally  we 
generate  a  frequency  comb  by  pumping  a  cavity  resonance  of  a  Si;,N4  microring  resonator  with  a 
single-frequency  tunable  diode  laser  amplified  by  a  ytterbium-doped  fiber  amplifier.  We  use  a 
polarizing  beam  cube  and  a  half-wave  plate  to  set  the  input  polarization  to  quasi-transverse 
electric  (TE)  and  a  bandpass  filter  to  reduce  the  amplified  spontaneous  emission  (ASE)  noise  of 
the  amplifier.  In  order  to  couple  light  into  the  chip,  we  focus  the  beam  onto  the  nanotaper 
coupler  of  the  S13N4  coupling  waveguide.  The  light  then  couples  into  a  ring  resonator  with  cross 
section  of  690x1400  nm,  coupling  gap  of  380  nm,  and  ring  radius  of  46  pm.  We  collect  the 
output  with  lensed  fiber  and  record  the  frequency  comb  spectrum  on  an  optical  spectrum 
analyzer  while  monitoring  10%  of  the  output  power  on  a  detector.  To  increase  the  circulating 
power  in  the  resonator,  we  tune  the  pump  wavelength  into  resonance,  which  causes  parametric 
oscillation  and  cascaded  FWM  to  occur,  resulting  in  comb  generation. 

We  experimentally  demonstrate  a  broadband  frequency  comb  spanning  830  to  1540  nm, 
shown  in  Fig.  17(b),(c),  which  represents,  to  our  knowledge,  the  broadest  comb  demonstrated 
near  the  visible  wavelength  range  [67].  This  wavelength  range  includes  both  the  D1  and  D2 
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atomic  absorption  lines  of  cesium  at  895  nm  and  852  nm,  respectively,  to  which  the  frequency 
comb  can  be  locked.  This  frequency  comb  also  encompasses  much  of  the  range  desired  for  OCT 
and  other  biological  imaging  applications.  The  demonstrated  comb  is  generated  with  an  amplifier 
power  of  10  W,  but  this  power  can  be  drastically  reduced  with  improved  coupling  efficiency  and 
quality  factor  (0.  For  the  device  demonstrated,  we  measure  a  coupling  efficiency  of  6dB  into 
the  chip  and  loaded  Q  of  500,000,  both  of  which  can  be  further  improved.  The  Q  is  lower  than 
expected  due  to  resist  polymerization  issues  during  fabrication.  Since  the  optical  parametric 
oscillation  threshold  is  inversely  proportional  to  the  square  of  the  Q,  an  increase  in  Q  will 
dramatically  reduce  the  required  amplifier  power.  These  improvements  will  lower  the  required 
amplifier  power  and  increase  the  bandwidth  of  the  comb.  Although  we  have  not  yet  seen 
evidence  of  single-pulse  modelocking  behavior,  simulations  confirm  that  this  can  occur  at  higher 
powers  or  with  the  fabrication  improvements  mentioned  above. 
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Figure  17.  (a)  Frequency  comb  simulation  for  1  pm  pump,  (b)  Experimentally  measured  comb 
spectrum,  and  (c)  comb  lines  generated  at  the  short  wavelength  edge  of  the  comb  indicated  in  inset 
of  (b). 


Next,  we  circumvent  the  traditional  challenges  in  achieving  microcombs  in  the  visible  range 
and  instead  simultaneously  employ  the  second-  and  third-order  nonlinearity  of  silicon  nitride 
(SiaNzO  to  generate  frequency  comb  lines  at  visible  wavelengths.  Silicon-based  materials  have 
traditionally  been  limited  to  third-order  (/l3))  nonlinear  interactions  as  they  are  inherently 
centrosymmetric  in  either  crystalline  (c-Si)  or  amorphous  (a-Si  or  SiNx)  forms.  Second-order 
nonlinear  optical  processes  require  breaking  the  symmetry  of  the  material.  The  recent 
demonstrations  of  S13N4  as  a  second-order  nonlinear  (/2))  material  [68,69,70,71,72]  has 
expanded  the  range  of  possible  applications  for  this  CMOS-compatible  platform.  The  symmetry¬ 
breaking  in  these  otherwise-centrosymmetric  materials  has  been  attributed  to  1)  high  levels  of 
strain  in  the  waveguide  material  [73],  2)  surface  effects  from  waveguide-cladding  interfaces, 
which  removes  the  bulk  symmetry  of  the  amorphous  film  [68,70],  and  3)  embedded  silicon 
defects  or  nanoclusters  [71,72],  The/(2)  coefficient  has  been  estimated  to  be  as  high  as  5.9  pm/V 
in  SiaN4  [72].  Second-order  nonlinear  processes  such  as  second-harmonic  or  sum-frequency 
generation  inherently  involve  wide,  octave-spanning  frequency  conversion  away  from  the  pump 
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frequency,  as  opposed  to  third-order  processes  such  as  FWM,  which  typically  involve  narrow 
conversion  near  the  pump  frequency.  Here  we  use  this  wide-band  attribute  of  /(2)  nonlinearity  to 
achieve  wide  frequency  conversion  from  near-IR  into  visible  wavelengths. 

We  generate  a  parametric  frequency  comb  in  the  near-IR  using  /3)  and  frequency-  translate  it 
into  the  visible  range  [Fig.  18(a)]  using  '/2)  in  an  integrated  ring  resonator  cavity.  A  single¬ 
frequency,  continuous- wave  (CW)  pump  laser  in  the  near-IR  is  coupled  into  a  high- Q  Si3N4  ring 
resonator,  which  acts  to  enhance  the  pump  power  and  serves  as  the  nonlinear  medium.  The  pump 
laser  first  undergoes  /(2)  second-harmonic  generation  (SHG),  producing  one  line  in  the  visible 
range.  As  the  cavity  circulating  power  increases,  the  pump  undergoes  /3)  FWM  optical 
parametric  oscillation,  which  produces  a  frequency  comb  in  the  near-IR  range  [Fig.  19(b)].  The 
comb  lines  then  mix  with  the  strong  near-IR  pump  through  a  y2)  sum-frequency  generation 
(SFG)  process  [Fig.  18(c)],  which  produces  a  set  of  comb  lines  in  the  visible  range  that  are 
generated  simultaneously  with  the  near-IR  comb.  The  generated  comb  lines  therefore 
simultaneously  address  two  widely  different  wavelength  bands.  Additionally,  all  processes  are 
achieved  in  a  single  integrated  device  using  a  single  CW  laser  source. 
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Figure  18.  (a)  A  single  near-IR  pump  is  coupled  into  a  microring  cavity  and  undergoes 
second-harmonic  generation,  near-IR  frequency  comb  generation,  and  sum-frequency  generation 
to  generate  simultaneous  near-IR  and  visible  wavelength  frequency  comb  lines,  (b)  Nonlinear 
photon  interaction  energy  conservation  for  degenerate  FWM  optical  parametric  oscillation  in  the 
near-IR.  (c)  Nonlinear  photon  interaction  energy  conservation  for  sum-frequency  generation 
between  one  pump  photon  and  one  comb  line  photon  generating  a  comb  line  in  the  visible. 


We  engineer  the  waveguide  dispersion  to  simultaneously  phase-match  near-IR  frequency 
comb  generation,  second-harmonic  generation,  and  sum-frequency  generation.  The  high  finesse 
of  the  Si3N4  cavity  enhances  the  efficiency  of  nonlinear  processes.  The  use  of  high-index 
contrast  waveguides  provides  the  ability  to  readily  control  the  waveguide  dispersion,  and  using 
Si3N4  eliminates  the  negative  effects  of  free  carriers  and  two  photon  absorption  at  the  operating 
wavelengths  [74].  Through  dispersion  engineering,  we  simultaneously  satisfy  the  phase¬ 
matching  conditions  for  each  nonlinear  process.  For  efficient  SHG,  the  effective  indices  of  the 
fundamental  and  second  harmonic  (SH)  modes  must  match.  This  phase-matching  is  achieved  by 
using  a  higher-order  waveguide  mode  for  the  SH  light,  as  in  [68].  The  third-order  TE  mode  at 
twice  the  fundamental  frequency  is  designed  to  have  the  same  effective  mode  index  as  the 
fundamental  TE  mode  of  the  pump.  We  use  a  finite-element  mode  solver  to  calculate  effective 
indices  of  each  mode  by  taking  into  consideration  material  dispersion,  a  4°  sidewall  angle 
introduced  during  fabrication,  and  the  waveguide  bending  radius  of  104  pm.  As  shown  in  Fig. 
19(a),  this  phase-matching  condition  is  satisfied  at  a  fundamental  wavelength  of  1540  nm  for  a 
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waveguide  with  dimensions  700x1400  nm.  This  waveguide  dimension  is  also  chosen  to  optimize 
optical  parametric  oscillation.  The  phase-matching  condition  for  optical  parametric  oscillation  is 
largely  dependent  on  group-velocity  dispersion  (GVD).  The  optimal  condition  for  oscillation  is  a 
low  level  of  anomalous  GVD  at  the  pump  frequency  [8].  This  dispersion  is  then  compensated  by 
the  nonlinear  phase  shift  in  order  to  achieve  the  proper  phase-matching  for  degenerate  FWM. 
As  shown  in  Fig.  19(b),  the  waveguide  dimensions  that  satisfy  SHG  phase-matching  also  exhibit 
a  sufficiently  low  level  of  anomalous  GVD  for  efficient  parametric  oscillation. 


Figure  19.  (a)  The  effective  mode  index  dispersion  of  the  fundamental  TE  mode  in  near  IR 
spectral  range  and  the  third  order  TE  mode  in  the  visible  spectral  range.  One  can  see  that  around 
the  pump  wavelength  of  1540  nm  (red  X),  phase-matching  occurs,  (b)  The  GVD  of  the 
fundamental  TE  mode  shown  in  (a).  The  low  level  of  anomalous  GVD  at  the  pump  wavelength  of 
1540  nm  allows  for  efficient  frequency  comb  generation,  (c)  Conversion  efficiency  vs.  wavelength 
for  %(2)  SHG  and  SFG  processes,  estimated  using  Eq.  (12).  Phase-matching  bandwidth  for  SFG  is 
wider  than  for  SHG,  allowing  for  multiple  comb  lines  to  be  converted  into  the  visible  range. 

Once  phase-matching  is  ensured  for  both  FWM  and  SHG  processes,  the  final  SFG  process  is 
automatically  phase-matched,  resulting  in  simultaneous  phase-matching  for  all  desired  processes. 
The  phase-matching  of  the  SFG  process  involves  one  photon  from  a  comb  line  mixing  with  one 
photon  from  the  strong  pump.  Fig.  19(c)  shows  normalized  conversion  efficiency  versus 
wavelength  for  SFG,  which  follows  the  same  phase  matching  condition  as  SHG  only  with  two 
photons  allocated  for  the  pump.  The  phase-matching  requirement  for  SFG  is  less  stringent  than 
that  of  SHG,  which  provides  the  opportunity  to  convert  a  wider  bandwidth  of  near-IR  comb  lines 
into  the  visible  range.  This  bandwidth  is  determined  primarily  by  the  slope  of  the  phase- 
mismatch  A k  with  respect  to  wavelength,  which  in  turn  is  determined  by  the  group  index 
mismatch  of  the  fundamental  and  higher-order  modes.  The  choice  of  a  different  higher-order 
mode  with  a  group  index  closer  to  that  of  the  fundamental  mode  could  further  increase  the 
conversion  bandwidth.  Overall,  these  three  nonlinear  processes  are  simultaneously  phase- 
matched  and  are  able  to  occur  together  in  the  same  cavity  structure. 
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Figure  20.  (a)  Near-IR  cw  pump  laser  and  frequency  comb  generation.  Several  different  states 
of  the  frequency  comb  are  shown,  (b)  Visible  second-harmonic  generation  and  frequency  comb 
lines  for  corresponding  spectra  in  (a).  Inset  micrograph  shows  the  visible  light  generated  by  the 
device. 


We  measure  up  to  17  comb  lines  frequency-translated  from  near-IR  into  visible  wavelengths 
extending  from  765-775  nm.  We  use  a  Si;,N4  ring  resonator  with  radius  104  pm,  waveguide 
dimensions  700  x  1400  nm  and  a  coupling  gap  of  560  nm.  The  resonator  has  a  Q  =  1.5><1 06.  All 
nonlinear  processes  occur  inside  the  resonator  cavity,  and  are  enhanced  by  this  high  Q.  The 
pump  is  produced  by  a  single-frequency  tunable  diode  laser  amplified  with  an  erbium-doped 
fiber  amplifier  (EDFA).  The  amplified  spontaneous  emission  (ASE)  noise  from  the  EDFA  is 
filtered  using  a  1-nm  tunable  bandpass  filter,  and  a  polarization  controller  is  used  to  couple  to  the 
fundamental  quasi-TE  mode  of  the  waveguide  using  a  single-mode  lensed  fiber.  The  pump  is 
tuned  into  a  cavity  resonance  near  1540  nm  (Fig.  20(a),  black  curve).  At  low  circulating  power  in 
the  cavity,  the  second-harmonic  is  generated  at  the  expected  wavelength  of  770.3  nm,  as  shown 
in  the  black  curve  of  Fig.  20(b).  The  maximum  second  harmonic  conversion  efficiency  we 
measure  is  -29  dB,  in  which  ~1  mW  of  the  second  harmonic  has  been  generated.  As  the  pump  is 
tuned  into  the  resonance,  the  circulating  power  in  the  cavity  builds  up,  which  leads  to  the  onset 
of  near-IR  frequency  comb  generation  as  well  as  frequency  conversion  of  the  comb  lines  into  the 
visible  range.  A  micrograph  inset  shows  the  generated  visible  light.  Several  states  [75]  of  the 
generated  comb  are  shown,  with  comb  lines  generated  every  5,  3,  and  1  free-spectral  ranges 
(FSR’s).  In  the  filled-in  comb  state  (red),  we  observe  a  visible  comb  bandwidth  of  nearly  10  nm, 
covering  17  FSR’s.  There  is  equal  frequency  spacing  of  217  GHz  among  the  near-IR  comb  lines 
and  the  visible  comb  lines.  The  SFG  conversion  efficiency  is  approximately  -20  dB.  We  believe 
this  conversion  efficiency  is  higher  than  the  SHG  conversion  of  the  pump  due  to  a  lack  of  precise 
pump  phase-matching  for  this  resonator.  This  is  evident  in  the  blue  curve  in  Fig.  20(b)  where  we 
see  that  the  converted  comb  line  five  FSR’s  away  has  a  higher  power  than  the  second-harmonic 
of  the  pump. 
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9.  Frequency  Comb  Generation  in  the  Mid-Infrared 

Silicon,  due  to  its  CMOS  compatibility,  wide  transparency  window,  and  high  third  optical 
nonlinearity,  is  an  ideal  platform  for  on-chip  comb  generation  deep  into  the  mid-IR.  However,  its 
linear  and  nonlinear  losses  have  until  now  prevented  the  realization  of  a  silicon  microresonator- 
based  comb  source.  Silicon  has  a  wide  transparency  window  from  1.2  pm  to  past  8  pm  [76]  and 
a  large  third-order  optical  nonlinearity  ( ri2  =  10'14  cm2/W  at  2.5  pm  wavelength  [77]),  which 
makes  it  an  excellent  platform  for  mid-IR  nonlinear  optics.  In  etched  silicon  microresonators, 
quality  factors  have  been  limited  by  scattering  losses  due  to  roughness  in  the  waveguide 
sidewalls,  which  is  made  worse  by  the  high  index  contrast  between  waveguide  core  and 
cladding.  The  dominant  nonlinear  loss  in  silicon  in  the  2.2  -  3.3  pm  region  is  three  photon 
absorption  (3PA)  -  a  process  where  three  photons  are  simultaneously  absorbed  to  excite  an 
electron-hole  pair  [78].  The  number  of  photons  lost  directly  to  3PA  is  small  (dominated  by  linear 
waveguide  losses),  but  the  generated  free-carrier  (FC)  population  will  induce  significant  optical 
losses19.  Here  we  overcome  both  silicon’s  high  linear  and  nonlinear  losses  to  demonstrate  a 
silicon  microresonator  frequency  comb  source,  and  to  show  the  fabrication  techniques  necessary 
to  achieve  an  on-chip  integrated  microresonator  comb  source  in  the  mid-IR. 

We  overcome  silicon’s  linear  losses  using  a  novel  etchless  fabrication  process  based  on 
thermal  oxidation  to  achieve  a  high  quality  factor  of  590,000  in  a  silicon  microresonator  at  a 
wavelength  of  2.6  pm.  The  silicon  waveguide  is  defined  in  an  etchless  manner  using  thermal 
oxidation,  instead  of  dry  etching,  to  form  the  waveguide  core  [79,80].  We  characterize  the 
resonator  quality  factor  by  sweeping  a  narrow  linewidth  source  across  the  resonance,  and  the 
loaded  quality  factor  is  220,000  and  the  one  is  590,000,  which  corresponds  to  a  propagation  loss 
of  0.7  dB/cm.  In  order  to  enable  broadband  comb  generation,  we  engineer  the  etchless 
waveguide  geometry  to  have  anomalous  group  velocity  dispersion  from  2.2  to  3pm.  To  mitigate 
silicon’s  nonlinear  loss,  we  embed  the  silicon  microresonator  in  a  reverse  biased  positive- 
intrinsic-negative  (PIN)  doped  junction  to  sweep  out  carriers  generated  from  three-photon 
absorption.  Significant  FC’s  accumulate  when  pumping  a  passive  silicon  waveguide  with  a 
continuous  wave  (cw)  laser,  limited  only  by  the  natural  FC  lifetime  of  the  structure.  Here,  we 
counteract  the  carrier  generation  while  using  a  cw  pump  by  extracting  the  generated  carriers 
using  a  PIN  junction  operated  in  reverse  bias  [81].  The  PIN  junction  prevents  the  electrical 
injection  of  carriers  into  the  waveguide  while  allowing  generated  free  carriers  to  be  swept  out  - 
with  effective  FC  lifetimes  demonstrated  as  short  as  12  ps  [82]. 

We  demonstrate  the  generation  of  an  optical  frequency  comb  between  2.1  pm  and  3.5  pm 
under  reverse  bias  conditions  (Fig.  21).  We  use  as  our  source  an  Argos  Model  2400  CW  optical 
parametric  oscillator,  which  is  a  tunable  from  2500  to  3200  nm  with  a  100  kHz  linewidth.  We 
couple  1.2  W  of  optical  power  at  2.59  pm  into  a  lens,  which  focuses  the  light  into  a  silicon 
nanotaper  on  the  chip.  The  input  coupling  loss  is  roughly  9  dB,  leading  to  150  mW  of  optical 
power  coupled  into  the  bus  waveguide  on  the  chip.  By  tuning  into  resonance  from  the  short 
wavelength  side,  we  achieve  a  soft  thermal  lock,  where  cavity  heating  is  compensated  by 
diffusive  cooling  to  maintain  constant  power  in  the  resonator.  Once  on  resonance,  we  extract  2.7 
mA  of  current,  using  10  volts  of  reverse  bias  applied  to  the  PIN  junction.  This  relatively  high 
current  is  an  indication  that  three-photon  absorption  is  significant  at  these  power  levels.  We 
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collect  the  output  light  from  the  nanotaper  using  another  lens,  and  couple  this  light  to  a  Fourier 
transform  infrared  spectrometer  (FTIR).  The  generated  frequency  comb  has  a  uniform  line 
spacing  of  127  ±  2  GHz,  well  below  the  7.5  GHz  resolution  limit  of  the  FTIR.  This  corresponds 
well  with  the  expected  free  spectral  range  (FSR)  of  the  resonator. 


Figure  21.  Broadband  mid-infrared  frequency  comb  generation  from  2.1  qm  to  3.5  qm  in  the 
etchless  silicon  microresonator.  This  frequency  comb  is  generated  with  150  mW  of  optical  power 
in  the  bus  waveguide,  and  with  10  V  reverse  bias  applied  on  the  PIN  junction. 


We  investigate  the  role  of  FC  absorption  in  silicon  microresonator-based  frequency  combs 
and  discuss  the  implications  of  operating  in  the  telecom  and  MIR  wavelength  ranges.  The 
nonlinear  dynamics  of  microresonator-based  frequency  combs  are  complex,  and  simulation  of 
the  process  over  broad  bandwidths  is  intensive.  We  apply  a  model  based  on  the  LL  equation,  but 
include  effects  of  multi-photon  absorption  and  FC  effects,  such  that 
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where  E(t,T)  is  the  field  within  the  microresonator,  Em  is  the  input  field  at  frequency  coo,  t  and  x 
correspond  to  the  “slow”  and  “fast”  time  variables,  respectively,  fa  is  the  round-trip  time,  ic  is  the 
power  transmission  coefficient,  a  is  the  round-trip  loss,  do  is  the  cavity  detuning,  L  is  the  cavity 
length,  /?£  is  the  k- th  order  dispersion  coefficient,  y  is  the  nonlinear  parameter,  /?2pa,  /?3pa,  and 
/?4pa  are  the  two-,  three-,  and  four-photon  absorption  coefficients  of  silicon,  respectively,  Atn  is 
the  effective  area  of  the  waveguide  mode,  and  a  and  /j  are  the  FC  absorption  (FCA)  cross  section 
and  FC  dispersion  (FCD)  parameter,  respectively  [83,84].  The  FC  density  Nc  is  governed  by  the 
rate  equation: 
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with  carrier  generation  rates  determined  by  multiphoton  absorption  and  the  recombination  rate 
determined  by  the  effective  FC  lifetime  reff .  Equation  (16)  describes  the  evolution  of  free  carriers 
on  the  “fast”  time  scale  within  a  round  trip  of  the  microresonator.  The  accumulation  of  free 
carriers  over  successive  round  trips  is  included  in  the  boundary  condition  Nc(t,  -Tr/2)  =  Nc(t  +  At, 
+-fe/2). 

We  consider  an  oxide-clad,  etchless  silicon  device  similar  to  that  used  in  the  experiment 
pumped  with  1  W  at  2.4  pm,  where  three-photon  absorption  (3PA)  becomes  the  dominant 
nonlinear  absorption  process.  The  spectra  generated  in  the  absence  of  3PA  (blue),  presence  of 
3PA  without  FC  (green),  and  presence  of  3PA  and  FC  (red)  are  presented  in  Fig.  22.  We  observe 
once  again  that  oscillation  does  not  occur  when  all  the  nonlinear  loss  terms  are  included. 
Interestingly,  in  contrast  to  the  telecom  frequency  comb,  it  is  evident  that  the  limiting  mechanism 
in  this  case  is  the  FCA  and  not  the  3PA,  as  simulations  which  include  3PA  without  FC  effects 
(green)  are  similar  to  those  without  any  nonlinear  loss  terms  (blue).  Since  the  rate  of  FC 
generation  is  intensity  dependent,  one  approach  to  mitigate  FCA  is  simply  to  lower  the  pump 
power  and  reduce  the  on-chip  optical  intensity  while  remaining  above  the  threshold  power  for 
parametric  oscillation.  The  threshold  power  is  primarily  determined  by  the  linear  losses  and 
effective  nonlinearity  of  the  resonator  and  is  on  the  order  of  a  few  mW’s  for  our  devices.  We 
note  that  the  threshold  power  is  independent  of  the  FC  lifetime  in  the  regime  presented  here  since 
oscillation  is  achieved  at  lower  powers  where  FCA  is  not  significant.  By  reducing  the  pump 
power  to  10  mW,  a  frequency  comb  is  generated  at  2.4  pm,  despite  the  presence  of  3PA  and 
FCA.  However,  this  approach  limits  the  achievable  parametric  bandwidth  and  is  not  conducive 
to  generating  a  broadband  comb. 


Wavelength  (nm) 


Figure  22.  (a)  Waveguide  GVD  and  (b)  simulated  frequency  comb  spectra  for  etchless  silicon 
microresonator  engineered  for  the  MIR  range.  Spectra  are  generated  in  the  absence  of  3PA  (blue), 
presence  of  3PA  with  no  FC  (green),  and  presence  of  3PA  and  FC  (red)  using  1  W  of  pump  power 
at  2.4  pm. 
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10.  Supercontinuum  Generation  in  Si-Based  Nano  waveguides 

Coherent  supercontinuum  generation  (SCG)  has  been  utilized  as  a  phase-coherent  broadband 
source  for  biological  imaging  and  molecular  detection  techniques.  A  phase-coherent  optical 
spectrum  is  also  critical  for  applications  including  pulse  compression  [85],  frequency  metrology 
[86,87],  and  wavelength  division  multiplexing  [88].  Additionally,  coherent  supercontinuum  (SC) 
with  an  octave- spanning  bandwidth  is  highly  desirable  for  the  detection  of  the  carrier  envelope 
offset  frequency  (fceo )  of  a  modelocked  laser  through  a  self-referencing  scheme  using  f-2f 
interferometry  [86,87],  enabling  a  fully-stabilized  frequency  comb  source.  For  some 
applications,  such  as  dual-comb  spectroscopy,  where  the  resolution  is  determined  by  the  comb 
spacing,  the  desired  source  repetition  rate  is  determined  by  the  sample  that  is  characterized 
[89,90].  Thus  a  SC  platform  that  can  operate  over  a  wide  range  of  pulses  energies  and  repetition 
rates  is  highly  advantageous.  Chip-based  SCG  offers  the  potential  for  cost-efficient,  large-scale 
production,  and  more  importantly,  the  potential  for  a  completely  integrated  photonic  SC  source. 
Current  CMOS-compatible  chip-based  SCG  platforms  include  silicon  [91,92,93],  amorphous 
silicon  [94],  silica  [95],  high-index  doped  silica  [96],  silicon  nitride  [97,98,99,100,101].  Silicon- 
based  waveguides  in  particular  provide  material  compatibility  with  existing  complementary 
metal-oxide-semiconductor  (CMOS)  fabrication  technology  which  allows  for  wide  scale 
implementation  of  integrated  chip-scale  devices.  Previous  demonstrations  of  SCG  in  Si;,N4  with 
near-infrared  [97,98,99,100]  and  mid-infrared  [101]  pump  lasers  have  achieved  octave- spanning 
bandwidths,  but  the  coherence  of  the  generated  SC  spectra  has  not  been  characterized. 

We  characterize  the  spectral  coherence  by  interfering  the  the  SC  generated  from  a  single 
input  pulse  with  the  spectrum  generated  from  a  previous  pulse.  From  this  measurement  we  find 
that  the  generated  SC  is  highly  coherent  over  the  majority  of  its  bandwidth  and  is  suitable  for/-2/ 
interferometry.  We  achieve  coherent  broadband  SCG  in  a  SF,N4  ridge  waveguide  through  careful 
design  of  the  dispersion.  SCG  at  a  specific  pump  wavelength  is  accomplished  through 
modification  of  the  waveguide  cross  section  to  tailor  the  appropriate  contribution  from  the 
waveguide  dispersion.  The  region  of  anomalous  GVD  at  the  pump  wavelength  allows  for 
coherent  SCG  through  soliton  fission  when  pumping  with  adequately  short  pulses,  and  the  two 
zero-GVD  points  allow  for  dispersive  wave  emission  on  either  side  of  the  pump  wavelength 
which  serves  to  extend  the  coherent  bandwidth  of  the  generated  SC  [17]. 


600 


800  1000  1200  1400  1600 

Wavelength  [nm] 


Figure  23.  (a)  Experimental  SCG  for  a  7.5-mm-long  waveguide  with  690x900  nm  cross- 
section  and  37  pJ  of  coupled  pulse  energy,  (b)  Left  axis:  Individually  simulated  SC  spectra  and 
ensemble  average  (black)  for  such  a  waveguide.  Right  axis:  Simulated  coherence  of  the  generated 
spectra  is  unity  over  the  entire  bandwidth. 
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We  demonstrate  SCG  using  a  1-GHz  repetition  rate  SESAM-modelocked  diode-pumped 
Yb:CALGO  laser  [102,103].  This  laser  allows  for  a  relatively  easy  coherence  measurement  due 
to  the  short  delay  in  the  path  length  (30  cm)  between  adjacent  pulses.  The  laser  pulses  have  a 
hyperbolic  secant  intensity  profile  and  a  center  wavelength  of  1055  nm.  An  optical  Faraday 
isolator  is  used  after  the  laser  to  protect  against  back  reflections  from  the  waveguide  facet,  and  a 
grating  pair  is  used  to  compensate  for  the  GVD  introduced  by  the  isolator.  Residual  GVD  and 
higher-order  dispersion  result  in  stretching  of  the  laser  pulses  from  63-fs  pulses  before  the 
isolator  to  91-fs  pulses  after  the  isolator  and  grating  pair.  Back  reflections  from  the  waveguide 
can  be  suppressed  in  subsequent  designs  through  use  of  an  angled  waveguide  facet,  eliminating 
the  need  for  an  isolator.  The  setup  for  SCG  is  the  same  as  described  above.  The  Si;,N4  waveguide 
is  7.5-mm-long  with  a  690x900  nm  cross  section.  Figure  23(a)  shows  the  SC  generated  with  236 
pj  pulses  incident  on  the  waveguide  with  an  8-dB  coupling  loss.  The  measured  SC  is  generated 
with  37  pJ  of  coupled  pulse  energy.  We  theoretically  investigate  the  coherence  of  the  SC 
spectrum  by  performing  numerical  SCG  simulations  using  the  split-step  Fourier  method  to  solve 
the  generalized  nonlinear  Schrodinger  equation  including  third-order  nonlinearity,  higher-order 
dispersion,  and  self-steepening.  To  characterize  the  spectral  coherence,  the  input  pulses  are 
seeded  with  quantum  shot  noise  [104],  and  the  first-order  mutual  coherence  function  for  the 
simulated  SC  is  calculated  based  on  [105].  The  simulated  SC  of  the  7.5-mm-long  waveguide 
using  91-fs  time-bandwidth  limited  pump  pulses  and  24.5  pJ  of  coupled  pulse  energy  is  shown  in 
Fig.  23(b)  (left  axis).  The  simulated  spectrum  corresponds  closely  with  our  experimental 
spectrum  and  displays  a  high  degree  of  coherence  over  the  entire  spectral  bandwidth,  Fig.  23(b) 
(right  axis).  In  order  to  experimentally  characterize  the  coherence  of  the  generated  SC  we 
perform  a  spectral  interference  measurement  [105,106]  in  which  the  output  from  the  waveguide 
is  sent  to  a  free-space  asymmetric  Michelson  interferometer.  One  arm  of  the  interferometer  is 
chosen  to  provide  a  delay  corresponding  to  the  pulse  period.  This  allows  for  interference 
between  adjacent  pulses,  creating  interference  fringes  on  the  optical  spectrum  analyzer  when  the 
SC  generated  by  subsequent  pulses  are  coherent.  A  single  sweep  of  the  OSA  records  the 
ensemble  average  of  >109  interference  events,  from  which  the  first-order  mutual  coherence 
^(Ajmay  be  calculated  according  to  [105], 


V(A)[/1(A)+/2(A)] 

2[/1(A)/2(A)]1/2 


(14) 


where  I\(X)  and  h(X)  are  the  measured  spectral  intensities  from  each  arm  of  the  interferometer, 
and  V(A)  denotes  the  fringe  visibility,  which  is  given  by  the  maximum  and  minimum  fringe 
intensity  as  V(X)  =  [/max(A)-/min(A)]/[/max(A)+/min(A)]  [107]. 
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Figure  24.  (a)  Left  axis:  Experimental  spectral  interference  measurement  (blue)  shows  high 
extinction  fringes  over  most  of  the  SC  bandwidth.  Right  axis:  visibility  (black)  is  close  to  unity 
over  much  of  the  wavelength  range,  (b)  Magnified  view  of  interference  fringes  and  extracted 
visibility  at  700  nm  and  (c)  1400  nm.  (d)  Extracted  coherence  (red)  is  close  to  unity  over  much  of 
the  wavelength  range. 

For  the  coherence  measurement  the  temporal  delay  between  interfering  pulses  is  selected  to 
give  adequate  resolution  (fringe  spacing)  over  the  entire  wavelength  range.  Light  output  from  the 
interferometer  is  collimated  into  a  fiber  and  sent  to  the  OSA.  We  use  a  fiber  that  is  single-mode 
at  1550  nm  for  measurement  of  wavelengths  longer  than  1250  nm  and  a  fiber  that  is  single-mode 
at  1060  nm  for  measurements  performed  below  1250  nm.  This  smaller  core  fiber  ensures  the 
spatial  overlap  between  the  interfering  pulses  at  wavelengths  longer  than  the  cutoff  wavelength. 
The  spectral  interference  measurement  is  performed  in  100-nm  steps  with  optimization  of  the 
collimation  and  the  alignment  of  the  output  for  each  wavelength  region.  The  interference 
measurement  is  shown  on  the  left  axis  of  Fig.  24(a).  As  seen  in  Fig.  24(a)  (left  axis),  interference 
fringes  exist  over  the  entire  spectral  range  above  the  noise  floor.  The  right  axis  of  Fig.  24(a) 
gives  the  extracted  visibility,  which  is  near  unity  for  the  majority  of  the  SC  bandwidth.  For 
measurements  taken  below  the  cutoff  wavelength  of  the  1060  nm  single-mode  fiber  (~  890  nm), 
the  multi-mode  nature  of  the  fiber  leads  to  a  sensitivity  to  the  fiber  collimator  input  coupling. 
This  effect  is  seen  in  the  reduction  of  the  fringe  extinction  from  740  -  840  nm  and  in  a  decrease 
of  the  visibility,  especially  at  the  edges  of  the  measurement  range.  Over  the  wavelength  range 
from  840  -  1750  nm  the  average  visibility  is  0.94.  This  high  average  visibility  results  from  the 
high  degree  of  phase  coherence  that  exists  between  interfering  pulses.  As  seen  in  Eq.  (14),  the 
visibility  is  equal  to  the  first-order  mutual  coherence  with  the  inclusion  of  a  normalization  factor 
that  takes  into  account  any  differences  in  the  intensities  of  the  individual  interferometer  arms. 
Importantly,  we  measure  visibility  near  unity  at  spectral  regions  that  are  an  octave  apart,  which 
enables  f-2f  interferometry  and  detection  of  the/ceo  of  the  modelocked  laser  source.  A  magnified 
view  of  these  spectral  regions  at  700  nm  and  1400  nm  is  shown  in  the  insets  of  Fig.  24(b)  and 
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(c).  We  measure  an  average  visibility  of  0.9  over  the  wavelength  range  from  675  nm  -  725  nm, 
and  an  average  visibility  of  0.99  over  the  wavelength  range  from  1350  -  1450  nm.  For  some 
wavelength  regions  we  are  limited  by  the  signal-to-noise  ratio  of  the  measured  SC  spectrum  and 
thus  could  not  extract  the  visibility.  Lastly,  we  calculate  the  first-order  mutual  coherence.  The 
calculated  coherence  shown  in  Fig.  24(d)  closely  resembles  the  extracted  visibility  with  slight 
differences  due  to  the  normalization  factor.  In  the  spectral  region  surrounding  the  pump 
wavelength  we  observe  modulations  in  the  measured  intensity  in  the  delay  arm  and  reference 
arm  of  the  interferometer,  which  distorts  the  coherence  calculation. 
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Figure  25.  (a)  Supercontinuum  spectrum  obtained  with  37  mW  of  coupled  average  power  at 
1.025  GHz.  (b)  Microwave  spectrum  showing  the  pulse  repetition  rate /rep  and  the  CEO  beat 
frequencies /ceo,i  and/CEo,2  at  0.024  GHz  and  1.001  GHz,  respectively. 


Finally,  we  detect  the  frequency  comb  offset,  given  by  the  carrier  envelope  offset  (CEO) 
frequency,  of  the  SESAM-modelocked  diode-pumped  Yb:CALGO  laser  [108].  In  order  to  detect 
the  CEO  frequency  of  our  lasers,  the  waveguide  output  is  collimated  and  sent  into  a  quasi¬ 
common  path /-to-2/ interferometer.  The  spectral  portion  around  1360  nm  is  frequency  doubled 
in  a  periodically  poled  lithium  niobate  (PPLN)  crystal  and  then  collinearly  recombined  in  space 
and  time  with  the  680-nm  components  of  the  supercontinuum.  After  passing  through  a  bandpass 
filter,  the  beat  signals  are  detected  using  a  gigahertz  avalanche  photodiode  (APD)  and  displayed 
on  a  microwave  spectrum  analyzer  (MSA).  In  Fig.  25,  we  show  the  detected /ceo  beat  signals 
and  the  corresponding  SC  spectra  that  were  used  for  /-to-2/  interferometry.  The  depicted  data 
was  obtained  with  30  pJ  (36  pJ)  of  coupled  pulse  energy,  corresponding  to  337  W  (345  W)  of 
peak  power  for  the  100-MHz  (1-GHz)  laser.  A  signal-to-noise  ratio  (SNR)  of  >  30  dB  with  a 
FWHM  of  ~  2  MHz  is  obtained  for  the  CEO  beat  of  the  1-GHz  laser.  The  SNR  of  the  CEO  beats 
is  sufficient  for  additional  microwave  signal  processing,  i.e.  it  will  allow  us  to  perform  further 
characterization  of  the  CEO  noise  properties  in  the  near  future.  Furthermore,  the  fact  that  the 
CEO  beat  center  frequency  can  be  shifted  over  a  substantial  range  when  adjusting  the  laser  pump 
power  is  an  essential  requirement  for  ultimately  stabilizing  the  CEO  frequency  to  an  external 
radio-frequency  reference  by  providing  feedback  to  the  current  of  the  pump  diode. 
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